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Abstract 

Taking into account the recent progress in the understanding of the lepton 
flavour effects in leptogenesis, we investigate in detail the possibility that the 
CP- violation necessary for the generation of the baryon asymmetry of the Uni- 
verse is due exclusively to the Dirac and/or Majorana CP- violating phases in 
the PMNS neutrino mixing matrix f/, and thus is directly related to the low 
energy CP- violation in the lepton sector (e.g., in neutrino oscillations, etc.). We 
first derive the conditions of CP-invariance of the neutrino Yukawa couplings 
A in the see-saw Lagrangian, and of the complex orthogonal matrix R in the 
"orthogonal" parametrisation of A. We show, e.g. that under certain conditions 
i) real R and specific CP-conserving values of the Majorana and Dirac phases 
can imply CP-violation, and ii) purely imaginary R does not necessarily imply 
breaking of CP-symmetry. We study in detail the case of hierarchical heavy 
Majorana neutrino mass spectrum, presenting results for three possible types 
of light neutrino mass spectrum: i) normal hierarchical, ii) inverted hierarchical, 
and iii) quasi-degenerate. Results in the alternative case of quasi-degenerate in 
mass heavy Majorana neutrinos, are also derived. The minimal super symmetric 
extension of the Standard Theory with right-handed Majorana neutrinos and 
see-saw mechanism of neutrino mass generation is discussed as well. We illus- 
trate the possible correlations between the baryon asymmetry of the Universe 
and i) the rephasing invariant Jcp controlling the magnitude of CP-violation in 
neutrino oscillations, or ii) the effective Majorana mass in neutrinoless double 
beta decay, in the cases when the only source of CP-violation is respectively 
the Dirac or the Majorana phases in the neutrino mixing matrix. 



^ Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 Sofia, 
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1 Introduction 



Baryogenesis through Leptogenesis [T] is a simple mechanism to explain the observed baryon 
asymmetry of the Universe. A lepton asymmetry is dynamically generated and then con- 
verted into a baryon asymmetry due to {B + L)-violating sphaleron interactions [2] which 
exist within the Standard Model (SM). A simple scheme in which this mechanism can be 
implemented is the "seesaw" (type I) model of neutrino mass generation [3]. In its minimal 
version it includes the Standard Model (SM) plus two or three right-handed (RH) heavy 
Majorana neutrinos. Thermal leptogenesis [U El [6] can take place, for instance, in the case 
of hierarchical spectrum of the heavy RH Majorana neutrinos. The lightest of the RH Ma- 
jorana neutrinos is produced by thermal scattering after inflation. It subsequently decays 
out-of-equilibrium in a lepton number and CP-violating way, thus satisfying Sakharov's 
conditions [7]. In grand unified theories (GUT) the masses of the heavy RH Majorana 
neutrinos are typically by a few to several orders of magnitude smaller than the scale of 
unification of the electroweak and strong interactions, Mqut = 2 x 10^^ GeV. This range 
coincides with the range of values of the heavy Majorana neutrino masses, required for a 
successful thermal leptogenesis. 

Compelling evidence for existence non-zero neutrino masses and non-trivial neutrino 
mixing have been obtained during the last several years in the experiments studying os- 
cillations of solar, atmospheric, reactor and accelerator neutrinos [HI El [IDl [HI [12] • The 
currently existing data imply the presence of 3-neutrino mixing in the weak charged-lepton 
current (see, e.g., [T3]): 

3 

i^iL = ^ Uij Vji^, I = e,n, r, (1) 
i=i 

where are the flavour neutrino fields, Uji^ is the field of neutrino Uj having a mass rrij 
and U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix [H]. The existing 
data, including the data from the ^H /3-decay experiments [I5], show that the massive 
neutrinos Uj are significantly lighter than the charged leptons and quarks d. The see-saw 
mechanism of neutrino mass generation [3] provides a natural explanation of the smallness 
of neutrino masses: integrating out the heavy RH Majorana neutrinos generates a mass 
term of Majorana type for the left-handed flavour neutrinos, which is inversely proportional 
to the large mass of the RH ones. 

Establishing a connection between the low energy neutrino mixing parameters and high 
energy leptogenesis parameters has received much attention in recent years [HI [IHl [201 
[2ll[221[231[2ll[25l[2ni[27l[28l|29]. These studies showed, in particular, that the number 
of phenomenological parameters of the seesaw mechanism is significantly larger than the 
number of quantities measurable in the "low energy" neutrino experiments. 

In the present article we investigate the link between the high energy CP-violation 
responsible for the generation of baryon asymmetry through leptogenesis and the leptonic 

^The data from the '^H /3-dccay experiments [15] imply nij < 2.3 eV (95%C.L.). More stringent 
upper limit on uij follows from the constraints on the sum of neutrino masses obtained from cosmologi- 
cal/astrophysical observations, namely, from the CMB data of the WMAP experiment [TB] combined with 
data from large scale structure surveys (2dFGRS, SDSS) (see, e.g., [17]). 
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CP violation at low energy, which can manifest itself in non-zero CP- violating asymmetry in 
neutrino oscillations and, in an indirect way, in the effective Majorana mass in neutrino less 
double beta {{l3l3)ou-) decay, \{m)\ (see, e.g., [HUl IHT] ). 

It was concluded in a large number of studies of leptogenesis performed in the so-called 
"one-flavor approximation" , that no direct link exists between the leptogenesis CP- violating 
parameters and the CP-violating Dirac and Majorana phases in the PMNS neutrino mixing 
matrix, measurable at low energies. In particular, an observation of leptonic low energy CP- 
violating phases would not automatically imply a non- vanishing baryon asymmetry through 
leptogenesis in the one-flavour case. This conclusion, however, does not universally hold 
[321 [HHl [3l] . Moreover, as was shown in [35] (see also [36]) and will be discussed in detail in 
the present article, the low-energy Dirac and/or Majorana CP- violating phases in U, which 
enter into the expressions respectively of the leptonic CP-violating rephasing invariant 
Jcp [SZ], controlling the magnitude of the CP- violation effects in neutrino oscillations, 
and of the effective Majorana mass | (m) | [3HI ED EH] can be the CP- violating parameters 
responsible for the generation of the baryon asymmetry of the Universe. Consequently, the 
leptogenesis mechanism can be maximally connected to the low energy CP- violating phases 
in U: within the scenario under discussion, the observation of CP-violation in the lepton 
(neutrino) sector would generically ensure the existence of a baryon asymmetry. 

The possibility of a direct connection between the high-energy leptogenesis and the low- 
energy leptonic CP-violation is based on the fact that a new ingredient has been recently 
accounted for in the leptogenesis mechanism, namely, the lepton flavour effects [32 l [331131] . 
As we have indicated above, the dynamics of leptogenesis was usually addressed within 
the 'one-flavour' approximation. In the latter, the Boltzmann equations are written for the 
abundance of the lightest RH Majorana neutrino, A'^i, responsible for the out of equilibrium 
and CP-asymmetric decays, and for the total lepton charge asymmetry. However, this ap- 
proximation is rigourously correct only when the interactions mediated by charged lepton 
Yukawa couplings are out of equilibrium. Supposing that leptogenesis takes place at tem- 
peratures T ~ Ml, where Mi is the mass of A'^i, the 'one-flavour' approximation holds only 
for T ~ Ml > 10^2 Qgv. For Mi > lO^^ Qgy, i.e., at temperatures higher than lO^^ Q^y, 
all lepton flavours are indistinguishable - there is no notion of flavour. The lepton asym- 
metry generated in A'^i decays is effectively "stored" in one lepton flavour. However, for 
T ~ Ml ~ 10^^ GeV, the interactions mediated by the r-lepton Yukawa couplings come into 
equilibrium, followed by those mediated by the muon Yukawa couplings at T ~ Mi ~ 10^ 
GeV, and the notion of lepton flavour becomes physical. 

The impact of flavour in thermal leptogenesis has been recently investigated in detail in 
[101 E21 E31 El] and in [211 IH] including the quantum oscillations/correlations of the asym- 
metries in lepton flavour space [32] . The Boltzmann equations describing the asymmetries 
in flavour space have additional terms which can significantly affect the result for the final 
baryon asymmetry. The ultimate reason is that realistic leptogenesis is a dynamical pro- 
cess, involving the production and destruction of the heavy RH Majorana neutrinos, and 
of a lepton asymmetry that is distributed among distinguishable lepton flavours. Contrary 
to what is generically assumed in the one-single flavour approximation, the AL = 1 in- 
verse decay processes which wash-out the net lepton number are flavour dependent, that 
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is the lepton asymmetry carried by, say, electrons can be washed out only by the inverse 
decays involving the electron flavour. The asymmetries in each lepton flavour, are therefore 
washed out differently, and will appear with different weights in the flnal formula for the 
baryon asymmetry. This is physically inequivalent to the treatment of wash-out in the 
one-flavour approximation, where the flavours are taken indistinguishable, thus obtaining 
the unphysical result that, e.g., an asymmetry stored in the electron lepton charge may be 
washed out by inverse decays involving the muon or the tau charges. 

When flavour effects are accounted for, the flnal value of the baryon asymmetry is the 
sum of three contributions. Each term is given by the CP asymmetry in a given lepton 
flavour /, properly weighted by a wash-out factor induced by the same lepton number 
violating processes. The wash-out factors are also flavour dependent. In the present article 
we perform a detailed analysis of the indicated flavour effects in leptogenesis. We show 
that the low energy Dirac and/or MajoranaCP- violating phases in U can be responsible 
for the generation of the baryon asymmetry of the Universe. We study also in detail the 
possible correlations between the physical low energy observables which depend on the CP- 
violating phases in f/pMNS - the rephasing invariant Jcp and the effective Majorana mass 
in neutrinoless double beta decay and the baryon asymmetry. 

The paper is organized as follows. In Section 2 we briefly review the existing data 
on the neutrino mixing parameters and the phenomenology of the low energy Dirac and 
Majorana CP-violation in the lepton sector. We note, in particular, that searching for 
CP- violation effects in u- oscillations is the only practical way to get information about the 
CP- violation due to the Dirac phase in the neutrino mixing matrix U (Dirac CP- violation), 
and that the only feasible experiments which at present have the potential of establishing 
the Majorana nature of light neutrinos Uj and of providing information on the Majorana 
CPV phases in U are the the neutrinoless double beta decay experiments searching for the 
process {A, Z) {A, Z + 2) + e~ + e~ . In Sections 3-7 we present a detailed discussion of the 
possible connection between the CP-violation at high energy in leptogenesis and the CP- 
violation in the lepton sector at low energies. We flrst derive the conditions of CP-invariance 
of the see-saw Lagrangian, i.e. of the neutrino Yukawa couplings, taking into account that 
the light and heavy neutrinos with deflnite mass are Majorana particles, and thus have 
deflnite CP-parities in the case of exact CP-symmetry (Section 3). We review briefly the 
arguments, based on the "one-flavour" approximation, leading to the conclusion that the 
connection between leptogenesis CP-violating parameters and the CP-violating phases in 
the PMNS matrix generically does not hold. We next discuss the conditions under which 
the lepton flavour effects in leptogenesis become important and present a brief summary 
of the results obtained recently on these effects, which are used in our analysis (Section 
4). In Sections 5, 6 and 7 we investigate the possibility that the CP- violation necessary 
for the generation of the baryon asymmetry of the Universe is due exclusively to the Dirac 
and/or Majorana CP- violating phases in the PMNS matrix, and thus is directly related to 
the low energy CP- violation in the lepton sector (e.g., in neutrino oscillations, etc.). The 
case of hierarchical heavy Majorana neutrino mass spectrum is studied in detail in Section 
5, where we present results for three possible types of light neutrino mass spectrum: i) 
normal hierarchical, ii) inverted hierarchical, and iii) quasi-degenerate. In Section 6 results 
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in the alternative case of quasi-degenerate in mass heavy Majorana neutrinos are derived, 
while in Section 7 we discuss how the results on leptogenesis in Section 4-6 will be 
modified in the minimal supersymmetric extension of the Standard Theory with right- 
handed Majorana neutrinos and see-saw mechanism of neutrino mass generation. Section 
8 represents Conclusions. 



Neutrino Mixing Parameters and Low Energy Dirac 
and Majorana CP- Violation 



We will use the standard parametrisation of the PMNS matrix: 

/ C12C13 S12C13 Si3e \ 

U= -S12C23 - ci2S23Sne'^ C12C23 - si2S23Si3e*'^ S23C13 diag(l,e*~,e*~ 

\ S12S23 - Ci2C23Si3e*'^ -C12S23 - Si2C23Si3e*'' C23C13 / 



(2) 



where c, 



[0,7r/2], S = [0,27r] is the Dirac CP- violating (CPV) 
phase and and 021 and a^i two Majorana CPV phases [121 1^ . We standardly identify 



cos Oij , Sij 



sin 9ij, 9ij 



case |Am^| 



m2 



> 0, where Ami, drives the solar neutrino oscillations. In this 

|Am^|, and 6q being the i/-mass 



"0 



|Am^2l' ^23 = Oa and 9u = Oq, 



squared difference and mixing angles responsible respectively for atmospheric and solar 
neutrino oscillations, while 6*13 is the CHOOZ angle [S]. The existing neutrino oscillation 
data, including the recent results of the MINOS experiment [15], allow us to determine 



Ami, 



Am^l, sin 612 and sin 2^23 with a relatively good precision and to obtain rather 

SZl HE]). The best fit values and the 95% C.L. 
and sin^ 29o-i read: 



stringent limits on sin ^13 (see, e.g. 
allowed ranges of AttIq , sin^ 612 



AmW 



(|Ami|)BF = 2.5 X 10-^ eV^ 
(Am|)BF = 8.0 X 10"^ eV^ 
fsin^ ^io)rf = 0.31, 



(sin^ 2e23)BF 



2.1 X 10-3 eV^ < lAm^l < 2.9 x 10"^ eV^ 
7.3 X 10-5 eV^ < Ami < 8-5 x IQ-^ eV^ 



0.26 < sin2^i2 < 0.36 
sin^ 2^23 > 0.90 . 



(3) 
(4) 
(5) 
(6) 



A combined S-u oscillation analysis of the global neutrino oscillation data gives [lUl W7\ 

sin^ ^13 < 0.027 (0.041) at 95% (99.73%) C.L. (7) 



The neutrino oscillation parameters discussed above can (and very likely will) be measured 
with much higher accuracy in the future (see, e.g., [13]). 

Depending on the sign of Am\ = Am|^ = Am^g which cannot be determined from 
presently existing data, the i/-mass spectrum can be of two types: 

- with normal ordering mi < m2 < m^, Am\ = Am'^^ > 0, and 

- with inverted ordering m^ < mi < m2, Am\ = Am^2 < 0; ''^i > 0; 

where we have employed the standardly used convention to define the two spectra. De- 
pending on the sign of Am\, sgn{Am\), and on the value of the lightest neutrino mass 
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(i.e., absolute neutrino mass scale), min(mj), the v-mass spectrum can be 

- Normal Hierarchical: mi<^ m2 -^m^, m2=(AmQ)2~ 0.009 eV, m3=|Am^|^~ 0.05 eV; 

- Inverted Hierarchical: m-^ <^ mi < m2, with mi 2 ^ |Ami|2 ~ 0.05 eV; 

- Quasi- Degenerate: mi = m2 = = m, m| S> |Am^|, m ^ 0.10 eV. 

Determining the nature of massive neutrinos, obtaining information on the type of u- 
mass spectrum, absolute z/-mass scale and on the status of CP-symmetry in the lepton 
sector are among the fundamental problems in the studies of neutrino mixing [13] . 

2.1 CP and T Violation in Neutrino Oscillations 

Searching for CP-violation effects in u- oscillations is the only practical way to get infor- 
mation about Dirac CP-violation in the lepton sector, associated with the phase 6 in U. A 
measure of CP- and T- violation is provided by the asymmetries [IHl HSl [50], [37] : 

(8) 

For S-u oscillations in vacuum, which respect the CPT-symmetry, one has [37] : 

^ 4M,r)^_^(e,r)^^^p^vac^ ^^'^ = (9) 

JcP = Im{t/eif/^2f/:2t^;i} = ^Sin2^i2sin2023cos2^i3sin^i3sin(5, (10) 

= si„ (^l) + sin (^ij + sin (^ij ^ (11) 

Thus, the magnitude of CP-violation effects in neutrino oscillations is controlled by the 
rephasing invariant associated with the Dirac phase S, Jcp [SZ] . The existence of Dirac CP- 
violation in the lepton sector would be established if, e.g., some of the vacuum oscillation 
asymmetries A^^p^r^^ etc. are proven experimentally to be nonzero. This would imply that 
<^cp 7^ 0, and, consequently, that sin 6^13 sin 5 7^ E|. One of the major goals of the future 
experimental studies of neutrino oscillations is the searches for CP- violation effects due to 
the Dirac phase in U (see, e.g., [T3l [52]). 



2.2 Majorana CP- Violating Phases and (/?/?)ojy-Decay 

As is well-known, the theories with see-saw mechanism of neutrino mass generation [3] of 
interest for our discussion, predict the massive neutrinos uj to be Majorana particles. We 
will assume in what follows that the fields i^jix) satisfy the Majorana condition: 

C{iy,f = uj, J = 1,2,3, (12) 

Let us note that the oscillations in matter, e.g., in the Earth, are neither CP- nor CPT- invariant [5T] 
as a consequence of the fact that the Earth matter is not charge-symmetric (it contains e~, p and n, but 
does not contain their antiparticles) . This complicates the studies of CP-violation due to the Dirac phase 
5 in z^-oscillations in matter (Earth) (see, e.g., [S2])- The matter effects in i/-oscillations in the Earth to a 
good precision are not T- violating [37j . however. In matter with constant density, e.g.. Earth mantle, one 
has [37]: a!^'^"^ — J^pF™^, J™p = Jcp Rcp, where the function Rqp docs not depend on 623 and S. 
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where C is the charge conjugation matrix: C ^7^C = — 7J, = —C, = C ^. 

Determining the nature of massive neutrinos is one of the most formidable and pressing 
problems in today's neutrino physics (see, e.g., [131 ES])- If are proven to be Majorana 
fermions, getting experimental information about the Majorana CPV phases in U , 021 
and asi, would be a remarkably difficult problem. The oscillations of flavour neutrinos, 
ui — > uii and ui uii, 1,1' = e,/i,r, are insensitive to the phases 021,31 [121 EI]- The 
Majorana phases of interest 021,31 can affect significantly the predictions for the rates of 
(LFV) decays /i— s>e + 7, r— s>/i + 7, etc. in a large class of supersymmetric theories with 
see-saw mechanism of z/-mass generation (see, e.g., [SI])- 

In the case of 3-z/ mixing under discussion there are, in principle, three independent CP 
violation rephasing invariants. The first is Jcp - the Dirac one, associated with the Dirac 
phase S, we have discussed in the preceding subsection. The existence of two additional 
invariants. Si and S2 is related to the two Majorana CP violation phases in U. The 
invariants 5*1 and 5*2 can be chosen as@ EHl IHT] : 

Si = Im {U:, Ur2} , ^2 = Im {U:, Urs} . (13) 

The rephasing invariants associated with the Majorana phases are not uniquely deter- 
mined. Instead of 5*1 defined above we could have chosen, e.g., S[ = lm{UeiU*2} 
or S'{ = Im [Ufj_iU*2} , while instead of S2 we could have used 5*2 = lm{Ue2U*-^}, or 
52 = Im |t/^2 f^^s}- The Majorana phases 021 and 031, or 021 and 032 = (031 — 021), 
can be expressed in terms of the rephasing invariants thus introduced [HI]: cos 021 = 
1 — {S[)^ / {\UeiUe2\'^) , etc. The expression for cos 021 in terms of 5*1 is somewhat more 
cumbersome (it involves also Jcp) and we will not give it here. Note that CP-violation due 
to the Majorana phase 0:21 requires that both Si = Im {UriU*2} 7^ and Re {UriU*2} 7^ 0. 
Similarly, 5*2 = Im{t/*2^r3} 7^ would imply violation of the CP-symmetry only if in 
addition Re {U*2Ur3} ^ 0. 

The only feasible experiments which at present have the potential of establishing the 
Majorana nature of light neutrinos uj and of providing information on the Majorana CPV 
phases in U are the experiments searching for the neutrinoless double beta ((/3/5)o^-) decay, 
{A,Z) {A,Z + 2)+e- + e- (see, e.g., [30l|53l [39] ) . The (/?/?)ov-decay effective Majorana 
mass, |(?Ti)| (see, e.g., [30|), which contains all the dependence of the (/3/3)ojy-decay amplitude 
on the neutrino mixing parameters, is given by the following expressions for the normal 
hierarchical (NH), inverted hierarchical (IH) and quasi-degenerate (QD) neutrino mass 
spectra (see, e.g., [39]): 



Ami sin^^i2e^°2^ + W Am^ sin^ ^ige^^^^'^-^'') 



rrii ^ m2 ^ (NH), 

(14) 



\{m)\ = y'lAmil |cos^^i2 + e*"2^ sin^^ial , 777,3 < mi < m2 (IH), (15) 
|(77i)| = mlcos^ 012 + 6'"^' sin^^ial , 7^1,2,3 = 77i > 0.10 eV (QD) . (16) 



'^Thc expressions for the invariants 6*1^2 we give and will use further correspond to the Majorana con- 
ditions (fT2|) for the fields of neutrinos ia,-, which do not contain phase factors, see, e.g., [31] . 
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Obviously, depends strongly on the Majorana CPV phase(s): the CP-conserving 

values of 021 = 0, ±7r [SZ] , for instance, determine the range of possible values of | (m) | in the 
cases of IH and QD spectrum, while the CP-conserving values of (031 —021) = 032 = 0, ±7r, 
can be important in the case of NH spectrum. As is well-known, in the case of CP-invariance 
the phase factors 

r/21 = e*°^^ = ±1 , r/31 = e*"«^ = ±1 , r/32 = e*°^^ = ±1 , (17) 

have a very simple physical interpretation [FTHHO] : rji^ is the relative CP-parity of Majorana 
neutrinos Ui and u^, r]ik = {Vi^^YVk^^ ^ Vi^kf ^ ^ being the CP-parity of i'i(k)- 

As can be shown, in the general case of arbitrary min(mj), depends on the two 

invariants Si and 5*2 [31]. In the chosen parametrisation of the PMNS matrix, Eq. (j2]), 
|(m)| depends also on Jqp- 

The (/3/3)oj/-decay experiments of the next generation, which are under preparation at 
present (see, e.g., [S3]), are aiming to probe the QD and IH ranges of \{m)\. If the {PP)ou- 
decay will be observed in these experiments, the measurement of the (/3/9)oj/-decay half-life 
might allow to obtain constraints on the Majorana phase 021 [311 [381 [58] . 



3 The CP-Invariance Constraints 

In the next Section, we will summarize the arguments leading to the conclusions that the 
leptogenesis CPV phases, responsible for the generation of the baryon asymmetry, can 
indeed be directly connected to the low energy CPV phases in U and, correspondingly, to 
CP violating phenomena, e.g., in neutrino physics. The starting point of our discussion is 
the Lagrangian of the Standard Model (SM) with the addition of three heavy right-handed 
Majorana neutrinos Ni {i = 1, 2, 3) with masses M3 > M2 > Mi > and Yukawa couplings 
It will be assumed (without loss of generahty) that the fields A'^- satisfy the Majorana 
condition: 

C(Njf = Nj, J = 1,2,3. (18) 

We will work in the basis in which the Yukawa couplings for the charged leptons are 
flavour- diagonal. In this basis the leptonic part of the Lagrangian of interest reads: 

/:i^P(x) = Cccix) + Cy{x) + CUx) , (19) 

where Ccc{x) and Cy{x) are the charged current (CC) weak interaction and Yukawa cou- 
pling Lagrangians, while C^{x) is the mass term of the heavy Majorana neutrinos Nf. 

Ccc = - -^/Z(a;)7ai^zL(x)iy"Hx)+h.c., (20) 

Cy{x) = \uN;^{x)H\x)iJiL{x) + hiH\x)'h.{x)^iL{x)+\i.c., (21) 
C^{x) = -]^M,N',{x)Ni{x). (22) 

Here ipiL and Ir denote respectively the left-handed (LH) lepton doublet and right-handed 
lepton singlet fields of flavour / = e,/i, r, tpj^^ = {un II), W"{x) is the H^^-boson field. 
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and H is the Higgs doublet field whose neutral component has a vacuum expectation value 
equal to f = 174 GeV. Obviously, Cy{x) + C^{x) includes all the necessary ingredients of 
the see-saw mechanism. At energies below the heavy Majorana neutrino mass scale Mi, the 
heavy Majorana neutrino fields are integrated out and after the breaking of the electroweak 
symmetry, a Majorana mass term for the LH flavour neutrinos is generated: 

= v'^X^ M-^ X = U*mU^ , (23) 

where M and m are diagonal matrices formed by the masses of Nj and Uj, M = 
Diag(Mi,M2,M3), m = Diag(mi, m2, ma), Mj > 0, > 0, and U is the PMNS ma- 
trix. The diagonalisation of the Majorana mass matrix mj,, Eq. fl2^ . leads to the relation 
([T]) between the LH flavour neutrino fields un and the fields Uj^ of neutrinos with definite 
mass and, correspondingly, to the appearance of the PMNS neutrino mixing matrix in the 
charged current weak interaction Lagrangian £cc{x)- 

In what follows we will often use the well-known "orthogonal parametrisation" of the 
matrix of neutrino Yukawa couplings [59] : 

A = - VMi?v^[/t , (24) 

V 

where R is, in general, a complex orthogonal matrix, R R^ = R^ R = 1. 

Before discussing leptogenesis and the violation of CP-symmetry associated with it, 
it proves useful to analyze the constraints which the requirement CP-invariance imposes 
on the Yukawa couplings A^;, on the PMNS matrix U and on the matrix R. If the CP- 
symmetry is unbroken, the Majorana fields Nj and Uk have definite CP-parities (see, e.g., 
[301157] ) rjj^^^ = ±i, //fc*"^ = and transform in the following way under the CP-symmetry 
transformation: 

UcpN,{x)Ulp = r/f CP 70 iV,(x'), r^f''^ = ^pf = ±^, (25) 
UcpM^)UIp = vfloi^kix'), r,f'' = ipl = ±i. (26) 

The sign factors and p^, and correspondingly, the CP-parities of the heavy and light 
Majorana neutrinos Nj and z/^, are determined by the properties of the corresponding RH 
and LH flavour neutrino Majorana mass matrices [SO]- They will be treated as free discrete 
parameters in what follows. 

Obviously, the mass term £^(x), Eq. 0221) . is CP-invariant. The requirement of CP- 
invariance of the Lagrangian Cy{x), Eq. (l2Tll . leads, as can be shown, to the following 
constraint on the neutrino Yukawa couplings: 

A*, = A,, (r/f CP)* V V""*, 3 = 1.2,?,, I = e, /i, r, (27) 

where r/' and rj^ are unphysical phase factors which appear in the CP-transformations of 
the LH lepton doublet and Higgs doublet fields, ipiL^x) and H{x), respectively. These phase 
factors do not affect any of the physical observables and, for convenience, we will set them 
to i and 1: rj'' = i, t]^ = 1. Thus, in what follows we will work with the constraint 0271) in 
the form 

A*, = A,, pf , pf = ±1 , J = 1, 2, 3, / = e, p, r, (28) 
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where we have used Eq. fl25|) . Note that if the CP-parity of a given heavy Majorana neutrino 
Nj is equal to (— i), i.e., if = —1, the elements Xji, I = e, p, r, of the matrix of neutrino 
Yukawa couplings A will be purely imaginary. 

It follows from Eqs. fl^^ and that if CP-invariance holds, the Majorana mass matrix 
of the LH flavour neutrinos generated by the see-saw mechanism is real (in the convention 
for the various unphysical phase factors employed by us): ruu* = rriy. This leads to the 
following CP-invariance constraint on the PMNS matrix U [30] : 

f/r, = f/;iP-,, J = l,2,3, / = e,/i,r. (29) 

In the parametrisation ([2]) we are using these conditions imply that the Dirac phase 5 = 
TTq, q = 0,1,2,..., and that the Majorana phases should satisfy: a2i = Tiq', aai = nq", 
g',g" = 0,1,2,.... 

Using Eqs. fpH|) . fpH|) and fl2I?l) it is easy to derive the constraints on the matrix R 
following (in the convention we are using) from the requirement of CP-invariance of the 
"high-energy" Lagrangian of interest ( fT9l) : 

R*, = R,kpfpl j,fc = l,2,3. (30) 

Obviously, this would be a condition of reality of the matrix R only if p^p^ = 1 for any 
j, k = 1,2, 3. However, we can also have pfp^ = —1 for some j and k and in that case the 
corresponding elements of R will be purely imaginary. 

The preceding results lead to the following (perhaps obvious) conclusions. 

i) If CP-invariance holds at "high" energy, i.e., if the neutrino Yukawa couplings satisfy (in 
the convention used) the constraints fl28l) . it will also hold at "low" energy in the lepton 
sector, i.e., the elements of the PMNS matrix will satisfy (in the same convention) Eq. fl2^ . 

ii) If the CP-symmetry is violated at "low" energy, i.e., if the PMNS matrix does not satisfy 
conditions (129|) and the CC weak interaction is not CP-invariant, it will also be violated 
at "high" energy, i.e., the neutrino Yukawa couplings will not satisfy Eq. ( 128|) and the 
Lagrangian jCy{x) will not be CP-invariant. Obviously, it is not possible to use the matrix 
R in order to render the Yukawa couplings CP-conserving. 

iii) If CP-invariance holds at "low" energy, i.e., if the CC weak interaction (1201) is CP- 
invariant and the PMNS matrix satisfy conditions ( l29l) . it can still be violated at "high" 
energy, i.e., the neutrino Yukawa couplings will not necessarily satisfy Eq. ( l28l) and the 
Lagrangian Cy{x) will not necessarily be CP-invariant. The CP- violation in this case can 
be due to the matrix R. 

As we will see, of interest for our further analysis is, in particular, the product 

Pjkml = Rjk Rjm Uif. Ulm , k ^ m. (31) 

If CP-invariance holds, we find from Eqs. ( l29l) and ( l30l) that Pjkmi is real: 

P;kna = Pjkml (Pf )' (PIY {P^rf = Pjkml , lm{P,,ml) = . (32) 

Let us consider next the case when conditions (j29l) are satisfied and Re{U*i^Urm) = 0, 
k < m, i.e., U*i^Urm is purely imaginary. This can be realised, as Eqs. ([2]) and ( l29l) show. 
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for 5 = 7rq, q = 0,1,2,..., and p^Pm = i-^-, if the relative CP-parity of the hght 
Majorana neutrinos z/^ and z/^ is equal to (—1), or, correspondingly, the Majorana phase 
a^k = 7r(2g' + 1), g' = 0, 1, 2, .... If Jcp = and both Re{U*^Ur2) = and ReiU^^U^s) = 0, 
CP-invariance holds in the lepton sector at low energies. In order for CP-invariance to 
hold at "high" energy, i.e., for Pjkmi to be real, the product RjkRjm, in the case under 
consideration, has also to be purely imaginary, Re{RjkRjm) = 0. Thus, purely imaginary 
U*jJJrm 7^ and purely real RjkRjm 7^ 0, i.e., Re(f/*^f/^m) = 0, Ivn^RjkRjm) = 0, in 
particular, imply violation of the CP-symmetry at "high" energy by the matrix R. 



4 Leptogenesis and the Baryon Asymmetry 

In this Section we will repeat briefly the arguments leading to the conclusion that the 
connection between leptogenesis CP-violating parameters and the CP-violating phases in 
the PMNS matrix generically does not hold. We will therefore first restrict our discussion 
to the "one-flavour" approximation. Then, we will review the recent results on leptogenesis 
which take into account the lepton flavour effects [32l [331 El]- In this case, the baryon 
asymmetry depends explicitly on the parameters of the lepton mixing matrix U. 



4.1 The One Flavour Case 

We assume that the heavy Majorana neutrinos Ni have a hierarchical mass spectrum, 
M2,3 3> Ml, so that studying the evolution of the number density of A^i suffices. The 
final amount of {B — L) asymmetry can be parametrized as Yb-l = nB-i/s, where s = 
27r^5f,,T^/45 is the entropy density and (7* counts the effective number of spin-degrees of 
freedom in thermal equilibrium {g^, = 217/2 in the SM with a single generation of right- 
handed neutrinos). After reprocessing by sphaleron transitions, the baryon asymmetry is 
related to the {B - L) asymmetry by [60] = (12/37) (Yb-l)- 
One defines the CP asymmetry generated by A^^i decays as 



_ Ez[r(iVi - Hli) - r(iVi ^ Hii)] _ i ^ im fM, 



m,)]_ 1 ^ Im [{\\%] (Mf\ 



^^[T{Ni^mi) + T{N,^mi)] Stt^ [AAt]^, '\Ml 

The wavefunction plus vertex contributions are included giving g{x) ~ — for x ^ 1. 
Notice, in particular, that ei denotes the CP asymmetry in the total lepton charge (i.e., 
trace over the lepton flavour index). 

Besides the CP parameter ei, the final baryon asymmetry depends on a single parameter. 



mj H{Mi) 

where H{Mi) denotes the value of the Hubble rate evaluated at a temperature T = Mi, 
~ 3 X 10~^ eV and 

_ _ (AAt)ii.;^ 

-1 = -j^ (35) 
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is proportional to the total decay rate of the heavy RH Majorana neutrino A^i. Notice 
again that is obtained by computing the total decay rate of the A^^i, that is by summing 
over all the flavours. 

The asymmetry in the total lepton charge is then given by 

YL^—ri{mi) , (36) 
9* 

where rj (mi) accounts for the washing out of the total lepton asymmetry due to inverse 
decays. 

The final baryon asymmetry in the "one flavour approximation" depends always upon 
the trace of the CP asymmetries over flavours, ei, times a function of the trace over flavours 
of the decay rate of the RH neutrino A'^i. This is an unexpected result because it suggests, 
in particular, that the lepton asymmetry, say, in the electron lepton number can be washed 
out by inverse decays involving the second and/or the third family (which erase only the 
muon and/or tau lepton charges). 

Let us turn now to the issue of the connection between the CP violating parameters 
in leptogenesis and the low energy CP- violating phases in the lepton sector, i.e., in the 
PMNS matrix U . The CP asymmetry ei can be written in terms of the diagonal light and 
heavy Majorana neutrino mass matrices introduced earlier, m = Diag(mi, m2, 771,3) and 
M = Diag(Mi, M2, M3), and the orthogonal complex matrix R: 



3Mi ^'^[Y.ipp^T^^P^^h^ip^^P^^P 



1/3 1 



|2 



3Mi Ep-^^^ 



Hp 



(37) 



1/3 1 



where we have summed over all lepton flavours / = e, /i, r. In particular, if the elements 
of the matrix R satisfy CP-invariance condition (!30l) . i.e., if -Rip, p = 1,2,3, is real or 
purely imaginary, the leptogenesis CP-asymmetry ei = 0. In the top-down approach, 
working in the basis in which the matrix of charged lepton Yukawa couplings and the RH 
neutrino Majorana mass matrix are diagonal, the matrix of neutrino Yukawa couplings can 
be written as A = V^Diag(Ai, A2, X3)Vl and the low energy leptonic phases can arise from 
the phases in Vl, i.e. in the left-handed (LH) sector, in Vr, i.e. in the RH sector, or from 
both Vl and Vr. However, from AA^ = Diag(A?, A^, Af) V/j = M^/^Rm R^M^/^v^, one 
sees that the phases of R are related to those of Vr. This means that in the "one-flavour" 
approximation a future possible observation of CP-violating phases at low energies in the 
neutrino sector does not imply the existence of a baryon asymmetry. Indeed, low energy 
CP-violating phases might stem entirely from the LH sector, i.e. V^, and hence be irrelevant 
for leptogenesis which would be driven by the phases in R, i.e. of the RH sector. 
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4.2 Taking the Flavour Effects into Account 



As we mentioned in the Introduction, the 'one-flavour' approximation rigorously holds only 
when the interactions mediated by the charged lepton Yukawas are out of equilibrium, that 
is at T ~ Ml ^ 10^^ GeV. In this regime, flavours are indistinguishable and one may 
indeed perform a rotation is flavour space to store all the asymmetry in a single flavour. 
At smaller temperatures, though, this operation is not possible. The r (/i) lepton doublet 
is a distinguishable mass eig enstate for T ~ Mi < 10^^ j^ig^) GeV. Flavoured Boltzmann 
equations may be written down for the asymmetry in each flavour A; = B/3 — Li 
[32l [33l . For the case of hierarchical RH neutrinos, they read 



ilD + 7S,AL=l) ( ^ - 1 ) , (3^ 



dz sH{Mi 



ei{-iD + 7S,AL=l) - 1 - + 7w,AL=l 



1 Ni J \ ^ J 1 1, 



(39) 



where z = Mi/T, T being the temperature, and there is no sum over / in the last term in the 
right-side of Eq. ( l39l) . In the above equation Y/Vi is the density of the lightest right-handed 
neutrino A''i with mass Mi. Y^^ are defined as Fa, = Fb/3 — 11,, where Yl, are the total 
lepton number densities for the flavours / = e, /i, r and Yb is the total baryon density. One 
has to solve the Boltzmann equations for Ya, instead of for the number densities Yi of the 
lepton doublets ii, since (A/ = B/3 — Li) is conserved by sphalerons and by the other SM 
interactions. Furthermore, all number densities Y, normalization to the entropy density s 
is understood and Y^^ and Y^'^ stand for the corresponding equilibrium number densities; 
7£) is the thermally averaged total decay rate of A'^i and 7s,al=i represents the rates for 
the AL = 1 scattering processes in the thermal bath. Notice, in particular, that 7s,al=i 
contributes to the asymmetry, as was recently pointed out in [31]. The corresponding 
flavour- dependent rates for wash-out processes involving the lepton flavour / are 7^5 (from 
inverse decays involving leptons £1) and 7w al=1) while AL = 2 scatterings may be safely 
neglected. The matrix A, which appears in the wash-out term connects the asymmetries 
in the lepton doublets to the asymmetries in the charges A^ by Yi = Aim ^a„ • The 
values of its elements depend on which interactions, in addition to the weak and strong 
sphalerons, are in thermal equilibrium at the temperatures where leptogenesis takes place. 
Below 10^ GeV in the SM, A is given by [H] 



/ -151/179 20/179 20/179 \ 
A = 25/358 -344/537 14/537 . (40) 
\ 25/358 14/537 -344/537 / 

In what regards the charged leptons, in the SM only the interaction mediated by the r 
Yukawa coupling is in equilibrium between 10^ and 10^^ GeV, and the lepton asymmetries 



^The equations of motion should include the terms accounting for the quantum oscillations among the 
different flavours [32j . They become relevant only in the proximity of the transition between the one-single 
flavour and the two-flavour state when the t flavour becomes distinguishable. 
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and S/3 — Li asymmetries in the e and /i flavour can be combined to Y2 = Ye+n and 



= ^Ae+A„- In this temperature range, A is given by 



A 



-920/589 120/589 
30/589 -390/589 



Above 10^^ GeV, all asymmetries can be combined to Yb-Li and A is given by A 
The asymmetry in each flavour is given by 



€-1 



3M: 



1 ^^(^ 



1/2 3/2 



e,/i,r. 



(41) 



-1. 



(42) 



1/3 1 



Obviously, the trace over the flavours of e; coincides with ei. It should be clear from 
Eq. ([32]) that we can have e; 7^ even if, e.g., i? is a real matrix and R ^ 1. If, however, 
i? is a diagonal matrix, e.g., if i? = 1, all three lepton number asymmetries vanish: ei = 0, 
I = e, jj,, T. 

Similarly, one has to deflne a "wash-out mass parameter" for each flavour / [321 [33] : 



nil 



3 X 10-3 eV 



\\ii\'v ' 
Ml 



nil 



T{Ni ^Hl) 



e,/i,r. 



(43) 



The quantity mi parametrizes the decay rate of A^i to the leptons of flavour /. The trace 
mi coincides with the irfi parameter deflned in the previous section. 

What is more relevant is that the total baryon asymmetry is the sum of each individ- 
ual lepton asymmetry. In the rest of the paper we will be concerned with temperatures 
(10^ ^ T ~ Ml ^ 10^^) GeV. In this range only the interactions mediated by the r Yukawa 
coupling are in equilibrium and the flnal baryon asymmetry is well approximated by 
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B 



37 



^2V 



417 
589 



m2 



390 

589 



(44) 



where €2 = eg e„, m2 



V (nil) 



mp 



and 



mi 



8.25 X 10-3 eV 



0.2 X 10-3 eV 

mi 



-1.16 



(45) 



Notice that the wash-out masses m2 and rrv in Eq. ( 14^ are multiplied by some numerical 
coefficients which account for the dynamics involving the lepton doublet asymmetries and 
the asymmetries stored in the charges A; = (1/3)5 — Li [H^ . 

From the generic expression for the baryon asymmetry, we deduce that the CP asym- 
metry in each flavour is weighted by the corresponding wash-out parameter. Therefore, the 
total baryon number is generically not proportional to ei. 
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The dependence of on the PMNS matrix elements is such that nonvanishing low energy 
leptonic CP-violating phases imply, in the context of leptogenesis and barring accidental 
cancellations, a nonvanishing baryon asymmetry. In the general case of complex matrix 
R, the low energy phases in the neutrino mixing matrix U could stem from the phases in 
the left-handed sector, in the RH sector, or in both sectors. This result only follows when 
the lepton flavour effects are correctly taken into account in the Boltzmann equations. In 
previous analyses of leptogenesis ignoring the flavour effects, the observation of low energy 
CP violation did not automatically imply the existence of a baryon asymmetry, since the 
possibility existed that the low energy phases could stem exclusively from the left-handed 
sector and hence be irrelevant for leptogenesis. 

This conclusion is by itself already quite interesting. However, we can go even further. 
Let us consider, for instance, the case in which the CP parities of the heavy and light 
Majorana neutrinos are such that pf = pj = 1 for all i,j = 1,2,3. In such a case. 



CP invariance corresponds to having all the elements of the matrix R real, see Eq. 
and 6 = 021 = 031 = (mod 27rjj. In the top-down approach, a real matrix R would 
correspond to the class of models where CP is an exact symmetry in the RH neutrino sector. 
The reason for this can be more easily understood working in the basis where the charged 
lepton Yukawa coupling and the RH mass matrix are diagonal, so that the neutrino Yukawa 
matrix is the only coupling in the leptonic Lagrangian that violates CP. More speciflcally, 
since the neutrino Yukawa coupling can be written in its singular value decomposition, 
A = V^Diag(Ai, A2, \3)Vl, CP violation in the RH neutrino sector is encoded in the phases in 
Vr, that can be extracted from diagonahzing the combination AA^ = V^Diag(A^, X^, A|)Vr. 
On the other hand, as we have seen, using the parametrization of the Yukawa coupling flMl) . 
this same combination of matrices can be written as AA^ = M^/"^ RmR'^ M^/"^ / . Comparing 
the two expressions it is apparent that R is real if and only if Vr is real, i.e. when there is 
no CP violation in the RH sector. It has been recently pointed out that the case of a real 
R matrix is naturally realized within the class of models based on sequential dominance 
[61] . In the case of R real, the flavour CP asymmetries and the baryon asymmetry depend 
exclusively on the phases of the left-handed sector, that are in turn uniquely determined by 
the low energy phases. Consequently, for real matrix i? 7^ 1, the leptogenesis mechanism 
is directly connected to the low energy CP-violating phases in f/pMNS- This connection is 
more apparent from the expression of the flavour CP asymmetries in the parametrization 

m 



€-1 



3Mi J2p ^/^pRlpUlp 

3Mi E/3 Ep>/3 ^rnpmp{mp - mp)RipRiplm. [UlpUip] 

12 



(46) 



If RijsRip is purely imaginary, RipRip = ±z |_Ri^i?ip|, we would get ±(mp -|- 



^ Alternatively, one may consider the case in which the CP parities of the heavy and light neutrinos 
are such that CP invariance is realized for the elements of the matrix R purely imaginary and 6^0, 
a2i = Q!3i = TT (mod 2tt). 
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mp)\RifsRip\Re{U*f^Uip) instead of (m^ — mp)RipRiplTci{U*pUip) in Eq. (H6|) . Purely imag- 
inary RipRip and real {UipUip) implies violation of CP-invariance. In order for the 
CP-symmetry to be broken at low energies, we should have both Re(?7f^f/;p) 7^ and 
l^UtpUip) ^ 0. 

In the next Section we will study in greater detail the possibility that the baryon 
asymmetry stems only from the low energy measurable CP violating phases in the neutrino 
mixing matrix U . 



5 Baryon Asymmetry from Low Energy CP- Violating 
Dirac and Majorana Phases in t/pMNS- Neutrinos 
with Hierarchical Mass Spectrum 

In what follows we shall assume that the matrix R has real and/or purely imaginary 
elements and that the heavy Majorana neutrinos possess a hierarchical mass spectrum, 
Ml < M2 < M3, with Ml having a value in the interval of interest, 10^ GeV ^ Mi < 10^^ 
GeV. We will investigate the case when the RG running of rrtj and of the parameters in U 
from Mz to Mi is relatively small and can be neglected. This possibility is realised (in the 
class of theories under discussion) for sufficiently small values of the lightest neutrino mass 
min(mj) [02], e.g., for min(mj) ^ 0.05 eV. The latter condition is fulfilled for the normal 
hierarchical (NH) and inverted hierarchical (IH) light neutrino mass spectrum. Under the 
indicated condition m^, and correspondingly Am\ and Amg, and U in Eqs. ( I46l) and 
Eqs. ( H3l) can be taken at the scale ~ Mz-, at which the neutrino mixing parameters are 
measured. 

Taking into account that in the case under discussion 62 = Ce + = — e,-, we can cast 
Eq. fl44p in a somewhat more convenient form: 

12 / /390_\ /417_\\ 
= - h ^"^2 , (47) 



37 V V589 J ' V589 J 
where m2 = rrie + rn^ and and ri{mi) are given in Eqs. (l43l) and ( l45l) . 



5.1 Normal Hierarchical Light Neutrino Mass Spectrum 

Given the inequalities mi -C m2 <^ m^, we will assume that the terms oc ^Jnii give sub- 
leading contributions to the lepton fiavour asymmetries e^, Eq. 0461) . and to the mass param- 
eters mi related to the wash-out effects, and we neglect them with respect to those oc ^Jfn^ 
giving the dominant contribution. This requires that ^mi|i?ii| <C ^^1x2 1 -R12 1 , \/'^?, \ R13 \ ■ 
In the indicated approximation, we find using Eq. 0461) and the fact that m2 = a/AttIq, 
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Im (e^^^^ U;,Uis) + 



Im {e'^''Re{U:,Ui3)) 



(48) 



SMiv^Ami /A 



■0 



Am^ 



Amg 
Am^ 



12 



+ I/? 



13 1 



xlm 



AuIq 



a/ Am 



A . 




13, 



,(49) 



where /323 = /5i2 + /^i3 = ai'g(-Ri2-Ri3)5 = arg(i?ij). The phase /523 parametrises the effect 
of CP- violation due to the matrix R in the asymmetry e;. In the case of CP-invariance, 
P23 = or 7r/20 depending on whether lm{Uj'2Ui-i) = or Re(?7;*2?7i3) = 0, respectively, and 
e; = 0, / = e,fi,T. 

From Eq. ([2]), it is straightforward to obtain ^(e^^^^f/fgf/zs) and Im(e*^2='Re(f/,*2f/z3)). 
The expressions, e.g., for Im(e*^^'''t/;*2^i3) read: 



Im (e*^^^ f/:2?7e3) = -S12C13S13 sin (^6 - + (^23 



-C23S23Cl2Sm(^— + (j23 



■S23'5i2Si3 sin 5 



Im {e'f^' f/;2f/.3) 



-Cl3 



C23S23C12 sm 



/«32 

V 2 



/^23 ) - C2,Si2Si3Sin (5 



V 2 

/tt32 

V 2 



f /523 
/323 



(50) 

)]' 

(51) 
(52) 



where 032 = a-si — a2i- Thus, e; depend on the same Majorana phase difference (or phase) 
the effective Majorana mass in (/5/5)oi/-decay, |(m)|, depends (see Eq. (fTSl) ) on. It follows 
from Eq. (H9l) that, as could be expected, for real or purely imaginary Ru and -R13 we have 



ee + + 



0. 



We are interested in the case in which the CP-violating phases in U play the role of 
leptogenesis CP-violating parameters. It follows from the preceding discussions that one 
has to consider two cases: i) P23 = (or more generally, (323 = ^rg, q = 0, 1,2, ...), and ii) 
P23 = 7r/2 (or more generally, P23 = (2g + l)vr/2, q = 0,1,2,...). We have to remember 
that if Im([/* ^«) 7^ (ReiU^^Uis) ^ 0) but ReiU^^Ui^) = (ImiUl^Uia) = 0), the case of 
/323 = , 71 {P23 = 7r/2, 37r/2) corresponds to violation of the CP-symmetry by the matrix R. 

Consider first the possibility of (323 = (vr). It is quite remarkable that, as it follows from 
Eqs. pOl) and (!52|) . for P23 = (vr) the asymmetry depends on the rephasing invariant S2, 



^More precisely, in the case of CP-invariance we have /?23 = ^rq or {2q + 1)tt/2, = 0, 1, 2, .. 
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Eq. f|T3|) : we have e,- oc Im {U*2Ur3) = 5*2. The requirement of a nonzero asymmetry, ^ 0, 
imphes S2 7^ 0, and correspondingly 032 7^ 2'7ik and/or {6 — 032/2) 7^ nk', k, k' = 0, 1, .... 
In order to reproduce the observed value of the baryon asymmetry, |sin(a32/2)| and/or 
I sin 6*13 sin((5 — 032/2)! have to be sufficiently large (see further). Since in the case under 
study P23 = or vr, the sign of er is not uniquely determined by the sign of Im {U*2Ur3)- 

In the alternative possibility of P23 = 7r/2 (37r/2) we have le,-! oc |Re(f/*2t/r3)|- Now 
^ provided 032 ^ vr(2A; + 1) and/or (6 - 032/2) ^ (7r/2)(2A;' + 1), k, k' = 0, 1, .... In 
this case | cos(a32/2)| and/or | cos^i3Cos(5 — 032/2)! have to be sufficiently large in order 
for leptogenesis to be successful. The sign of e,- is again not uniquely determined by the 
sign of Re {U*2Ur3). 

The maximal value of |Im(e*^^^f/*2^r3)| and of le,-!, is reached for /?23 = 0; n (/523 = 
7r/2; 3n/2) at | 032 = 7r(2A; + 1) and 5 = 2'Kk (032 = 2TTk and 5 = n{2k + 1)), A; = 0, 1, 

|lm {e'^^^ U;2Ur3) I = C23C13 (S23C12 + C23S12S13) ^ 0.47 , /523 = ^ g , g = 0, 1, 2, 



max 



2 

(53) 

where we have used the best fit values of sin^ 2^23 and sin^ given in Eq. ([5]) and the upper 
limit si3 < 0.2 (see Eq. ([ZD). For S13 = we get \lm{e'^^^U*2Ur3)\ ^ 0.42 | sin(a32/2)| < 0.42 
while if 032 = 0, one obtains \lm{e''^^^U*2Ur3)\ = 0.27 |si3sin5| < 0.054. Thus, the effect 
of the term oc |si3sin5| in 15*21, and correspondingly in |e,-|, can be significant only if 
I sin(a32/2)| < 0.20. Note that for (323 = 0; tt (/323 = vr/2; 37r/2) max (|Im(e^^2^f/;2?7r3)|) 
corresponds to Re{U^2Ui3) = (Im(f/f2t/i3) = 0), / = e,fi,T, i.e., to CP-conserving values of 
^32 and 6. Nevertheless, in this case we still have le,-! 7^ as a consequence of the violation 
of the CP-symmetry by the matrix R. 

We turn next to the dependence on the parameters in R. Taking into account that R 
is, in general, a complex orthogonal matrix, we get for the elements of R of interest, Rij, 
j = 1,2,3, obeying the CP-invariance constraints: 

(|i?n|Vr + 1^121^2 + 1^131^3) Pf = 1 • (54) 
The CP-asymmetry \€t-\ is proportional to the factor r: 

I-R12-R13I 



2 ~ 

^"^0 1 ^ I p 12 I I D 12 



(55) 



It is clear that as long as I-R12I ~ |-Ri3|; ^ will not act as a suppression factor for the 

asymmetry le^-l; the latter can be suppressed if, for instance, |i?i2| <C I-R13I or |i?i3| <^ 
{Aml/Aml)^Ri2\. 

We now focus on the wash-out effects. The mass parameters related to the wash-out 
effects for the three lepton number asymmetries read: 

2 

, / = e,/i,r. (56) 



mi ^ \/Aml 



2xi 



I R12 U12 + Ri3 



^Note that for /Jss = 7r/2; 37r/2 we have \lm{e^l^^^U*2Ur3)\ = \Rc{U*2Ur3)\, and the maximum of 
|Im(e''^=3J7*2C/r3)| coincides with the maximum of |Re(J7*2^i-3)|- 
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Using the unitarity of the PMNS matrix U we obtain: 



rrie + rriij, + 



AmQ \ 2 |2 I |2 



(57) 



Since mi > 0, I = e, /i, r, each of the individual wash-out mass parameters is hmited from 
above by the expression in the right-hand side of Eq. ( 157|) . 

Except for strong fine tuning, the contribution due to sin ^13 is subdominant and can 
be neglected at first order. Therefore, for simplicity, we take sin ^13 = in the following 
analysis and we get: 



7712 



rrir 



'Am\ 



/A 



■nrir 



( Am, 



I ® 

\^AmJ 



\R 



121 



2^2 

23 



1^12^13! S23C23Cl2COs((/3i3 -/?i2 



^32 . 

2 ■ 



a/A 



AttIq 

Am\ 

\ 
4 



I R |2^2 „2 
|-n:i2| C;^2'523 



I P |2^2 

C23 



'^(S^J 1^12^13! S23C23Cl2COs((/3i3 



^ a32. 



(5^ 



(59) 



It follows from Eq. (jHj) that the baryon asymmetry |Yb| can be zero if |?7(0.66r?v) — 
r;(0.71m2)| = 0, although je^l 7^ and |ee + 7^ 0. This corresponds to the physical case 
when the asymmetries generated in the lepton doublet charges r and (e + yu) are equal in 
magnitude, but have opposite signs. Such a possibility can occur if the following relation 
holds: 

'390_\ /417^ 



589 



rrtr 



V 



589 



m2 



One solution is given by: 



_ 390 _ _ 

1712 = mr = 0.935 mr . 

417 



(60) 



(61) 



For fixed I-R12P and I-R13P, Eqs. flFTI) and flUTl) determine the value of for which we can 
have \Yb\ = 0: 



x/A 



m. 



0.935 



Am|\ 2 

-T — 2" -^12 + -^131 
Am^ 



(62) 



Since is calculated from Eq. fl56p . given the neutrino oscillation parameters and I-R12I 
and I-R13I, the requirement that the value of r7v obtained from Eq. coincides with 
that determined by Eq. leads to a combined constraint on the quantities |^7r2p and 



2 Re(-Ri2f/*2-Ri3 ^ts) which depend on the phases 6 and 032. The phase factors associated 
with the latter vary only between (—1) and 1. Therefore it is not guaranteed that condition 
(16T|) can be satisfied and we can have \Yb\ =0 with le,-! 7^ and |ee-|-e^| 7^ for any possible 
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values of R12 and R13. Taking for simplicity — P12 — 032/2) = vr/2, a solution is given 
by: 



\R 



12 



0.94|i? 



, cos 29' 



23 



13 



1 - 2.069c?2si3 



(63) 



which holds for 2.069 c|2% ^ 1 and (cos2e23-0.069s^3)/(l-2.069c?2S^3) > 0. For example 
if c^2 = 0.69, it can be satisfied for ^ 0.48. 
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A general solution of Eq. fl60|) can also be found and is given by: 



589^ 



m2 ttIt- 



390 417 



lO"*' eV^ 



y-i 

,1.16 _ 



,0.08 



2 

2.16 



(64) 



where y = (417/390)m2/mT-. If, for instance, /3i3 — /3i2 — ^32/2 = 7r/2 and I-R13P dominates 
in and m2, y does not depend on the parameters in R and we can solve Eq. 



\R 



13 1 



10^6 eV^ 



589 589 / y 



,0.08 



ci3si3 390 417V|/i-i6_ r 



2 

2.26 



0.25 . 



(65) 



We will analyse next the dependence of the baryon asymmetry on the parameters in 
R. For simplicity, we take the values of 032 and of /3i2,/9i3 = 0,7r/2, which maximize the 
CP-asymmetry in Eq. (l49l) . The wash-out parameters read: 



m2 



Ami 



ArriQ 



12 



2 2 ' 

Cl2'^23, 



+ i^i3r4) 



Am; 



■0 



f |-Rl2|^C^2S23 + |-Rl3rC23 ) . 



2„2 



Am^ 



(66) 
(67) 



Case A: Strong wash-out. This case is realized if m2,r ^ 2 x 10 ^ eV. The latter 

2 o in-3„A//^ /Amlcj^sls) = 0.64, and/or for {Ri^l"^ > 



condition is satisfied for |i?i2|^ > 2 x 10 ^eV 

0.08. The baryon asymmetry can be approximated as: 
-4 ^\r. 1 ifi /Con o XX 1 n-4 



2 X 10-3 eV/(v/ Amicus) ~ 



\Y, 



B\ 



c- 



I-R12II-R13I 



f Am^ 
.Am\ 



2 X i 



\R 



12 



\R 



589 2 X 10-4 eVy ie /589 2 x 10"^ eV\i-i6 



390 



m,- 



V417 



m2 



. (6^ 
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The constant C is defined as 
9Mi 



C = 



Am\ 



Am; 



■0 



Am^ 



* C23S23C12 ^ 1.1 X 10 ^ 



Ml 



1011 GeV 



(69) 



where we have used the present best fit values of the neutrino oscillation parame- 
ters. If the l-RisP term in the wash-out factors dominates, i.e. if l-RioP/l-R 



1-13 



/ \ 1/2 

Ami/AmM ^23/(1 

~ '^12-^23) — "^-3; the dependence on the R parameters becomes: 



B\ 



2.7 X 10-3 C 





-R12 




-R13 


3.32 


„2.32 
''23 



390 c^3\i i6 
417 ^ 



(70) 
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For maximal atmospheric neutrino mixing there is a strong cancellation and the resulting 
baryon asymmetry has an additional suppression factor ~ 0.075: 

|lB|~4.5xlO-C^g|L. (71) 
Taking into account the dependence of on sin 

2^23, we find that for sin^ ^23 = 0.34 (0.66) 
\Yb\ is larger by a factor of 9 (11). Notice that the asymmetry changes sign when sin^6'23 
increases from 0.34 to 0.66. Thus, there is also a value of sin^ 6'23 in the interval [0.34,0.66], 
for which Yg = (up to corrections ~ S13). The asymmetry decreases very rapidly with 
|i?i3| and therefore the maximal asymmetry should correspond to relatively small l-Riaj. 

In the alternative case when I-R12P dominates in m^, which is realised for I-R12P ^ 
7.9 I-Risp, the asymmetry \Yb\ is proportional to: 

■Ami \ 1-08 



\Yb\ ~ 2.7 X 10-^ C ( =^ 
' ^ Ami 





-R13 




-R12 


3.32 



1-16 /417,, o 



'^12'^23 ) ^390^^ '^12'^23/ 



(72) 



In this case, \Yb\ is not suppressed for maximal atmospheric neutrino mixing. Varying the 
solar and atmospheric angles within their 95% CL. ranges induces a factor of a few in 
the predicted baryon asymmetry. Also in this case the baryon asymmetry is a decreasing 
function of |-Ri2|- 

Case B: weak wash-out. For sufficiently small values of I-R12I and I-R13I, one enters the 
weak wash-out regime. The asymmetry is given approximately by: 



I-R12II-R131 



/'Am?.-, \ 2 



390 m; 417 m2 



589 8.25 X 10-3 eV 589 8.25 x 10-3 eV 



(73) 



If the term oc I-R13P dominates in mT-2, we can have again a partial cancellation between 
the two leading terms in \Yb\ for maximal atmospheric neutrino mixing, ^23 = 7r/4: 

27 
390 



^4.0C |i?i2||i?i3| 



cos 26*23 - 7T7T7T''^23 



0.12C|i?i2||i?i3| . (74) 



The asymmetry \Yb\ is bigger by approximately a factor of 10 if the atmospheric mixing 
angle assumes the values sin^ ^23 = 0.34, 0.66. The asymmetry has opposite sign for these 
two values, indicating that \Yb\ goes through zero (up to corrections ~ S13) when sin^ ^23 
is varied from 0.34 to 0.66. 

If the term cx: I-R12P is the dominant one in r?V2, one finds: 

Am^ /417 27 \ 

|y,| . 4,0 C^lfldlftal (5^ + ^44) - O.UC\R,.A\Ru\ ■ (75) 

We notice that in the case of weak wash-out one finds that the asymmetry increases with 
the values of |i?i2| or I-R13I. 

On the basis of the above discussion, we can conclude that more than one local max- 
imum of \Yb\ is present and that the absolute one is found, in general, in a regime which 
interpolates between the weak and strong wash-out regimes. However, it is possible to show 
that in this region Yb can be even zero if Eq. (l60l) is satisfied. 



21 



5.2 The case of A^3 decoupling 



We will assume further for simplicity that |-RiiP <C II-R12P ± l-RisPl- This corresponds to 
the case of "decoupling" of the heaviest Majorana neutrino N3 [H31 EH [2H] • For P23 = 0; vr 
{P23 = 7r/2; 3n/2), which implies P2P3 = 1 ^6 get: 

\Ru\\l)\Rl3\' = P2 Pi ■ (76) 

Obviously, if P23 = 0; vr we have I-R12P + l-RisP = 1, while for P23 = 7i"/2; 37r/2, one finds 
|i?i2p — I-R13P = ±1. It is not difficult to find the maximal value of r for p2 Pi = 1: 

maxr = -(^^j =1.2, P^pf = l, /323 = -k , k = 0,1,2, ... (77) 



The maximum is reached when (323 = 0; vr at I-R12P = (1 + a/ ^ttlq/ a/ Am^) ^ = 0.85, 
while for P23 = 7r/2; 37r/2 it corresponds to |i?i2p = (1 - ^/ Ami / ^ Aml)-^ - In 
the latter case I-R12P — l-^isP = 1 ^ind I-R12I, I-Risl, in general, are not limited from above 
1^. However, it is not difficult to convince oneself that we can have successful leptogenesis 
only if |_Ri2| and I-R13I are not exceedingly large, namely for I-R12I, I-R13I ~ 10. Indeed, 
barring accidental cancellations, the wash-out mass parameters mi, I = e,fi,T, increase 
with the increasing of I-R12I (l-RisP = I-R12P 1); and for |i?i2| ~ I-R13I ~ 10, one would 
typically have m/ ~ 1 eV. The corresponding efficiency factors would be exceedingly small 
and \r] (f|§?^) — V (Hg^) I ~ 10~^ (extremely strong wash-out regime), which makes it 
impossible to reproduce the observed value of the baryon asymmetry for Mi ^ 10^^ GeV. 
The preceding rather qualitative analysis can obviously be refined to obtain more precise 
upper limits on |i?i2| and I-R13I satisfying I-R12P — l-RiaP = 1- 

If P23 = 7r/2; 37r/2 and P2Pi^ — ~1; one always has r < 1. In this case I-R12P < l-RisP- 
For I-R12I <^ I-Risl we have r -C 1 and the asymmetry \er \ will be suppressed by the r-factor, 
Eq. (l55l) . We will not analise this case further. 



Given the fact that v^Am^ = 0.05 eV, y^Am'^/^/Aml = 0.18, r < 1.2, Eq. ([771), and 
|Im(f/*2f/r3)| ^ 0.47, we obtain for the maximal value of the asymmetry |er| in the cases of 
interest {\Ru\^ (+)|i?i3p = 1, ^23 = {^/2)k, k = 0,l, 2, ...): 



, 3Mi J Ami , , s / a/ Ami \ / Mi \ 

< 0.19 (0.27) . 1.9 (2.7) x 10- j (^ j . (79) 

where the minus (plus) sign and the value 0.19 (0.27) correspond to P23 = tt/c, k = 0,1, 2, ... 
(/523 = (7r/2)(2A; + l), A; = 0,l,2,...). 



^Note that if I-R12P + ~ 1, one can obviously use the paranictrisation R12 = cosoj, R13 ~ sinoj, 

while in the case of |i?i2p — l^isP = 1; we have R12 — cos{iuj) = coshoj, i?i3 = sin(jcj) = isinhLj, uj being 
a real parameter. 
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Let us note that, although for the values of the parameters we have used the asymmetry 
le,-! has a maximal value, the baryon asymmetry \Yb\ does not necessarily has a maximum 
since the maximum of \er\ does not correspond, in general, to a maximum of the effective 
wash-out factor |?7(0.66m^) — 77(0.717712) | in the expression for \Yb\ (see Eq. ( l44l) ). As can 
be shown and Figs. [T] and [2] illustrate, however, for real R12 and -R13 {(323 = 0; it) satisfying 
I-R12P + l-RisP = 1) the maximum of |eT-| with respect to the parameter R12 practically 
coincides for S = {S = Tik, k = 0,1, ...) and CP-violation due only to the Majorana phase 
Q!32, with the maximum of IYbI- In the case of CP- violation generated only by the Dirac 
phase 6 (032 = 0), the maximum of |Yb| occurs for P23 = at slightly smaller value of I-R12I, 
namely at I-R12I — 0.86, compared to the value of I-R12I — 0.92 at which the maximum of 
le,-! takes place. We will work for convenience with R12 and -R13 maximising |eT-|, for which 
we have simple analytic expressions in terms of the ratio Am|/A?7i^. In most of the cases 
we will discuss these values of R12 and R13 maximise also \Yb\. 

Consider next the wash-out parameters mj, I = e, fi, r. In the case of real Ru and R13 
we have P23 = 0; vr and I-R12P + |-Ri3p = 1- Taking into account this constraint we get 
from Eq. [281E3]: 

^J AuIq <rrie + + rrLr < \J Am^ . (80) 

The maximum and the minimum are reached respectively for I-R12P = and I-R13P = 0, and 
in both cases |e^| = 0. If (323 = 7r/2; 3n/2, we have |-Ri2p — I-R13P = 1 and for me + m^^ + m(r 
we get the same lower bound as in Eq. flHUj) . but no upper bound, in general, since I-R12P 
is not limited from above. The requirement of successful leptogenesis leads to an upper 
bound roughly of the order of 1 eV. For I-R12P = (1 (-)\/^"^/ V ^"^a)~^! corresponding 
to the maximum of r and of |e^|, we find 



TUe + rrifj, + m-j- = 2 y Am| 

In what follows we will analyse the case of real R12 and R13. For real R12 and R13 we 
have P23 = 0; vr corresponding to sgn(i?i2-Ri3) = +1; (—1). 

5.2.1 Leptogenesis due to Majorana CP-Violation in ?7pmns 

We shall consider first the interesting possibility of CP-symmetry being violated by the 
Majorana phase 032 in U, and not by the Dirac phase 6. To be concrete, we choose 
5 = and real R12 and -R13 {(323 = 0; tt) which maximise r, \er\ and \Yb\, i.e., I-R12P = 
(1 + ^yAml/^yAml)-^ = 0.85, |i?i3p = 1 - |i?i2p = 0.15. For 032 = 0; 2tt, the CP- 
symmetry is not violated. If 0^32 takes the CP-conserving value of 032 = vr, CP-symmetry 
is violated by R and |e^| 7^ 0. It should also be noted that the terms oc sin 6*13 in the 
expressions for and the wash-out mass parameters mi are sub-dominant in the case 
being studied. 



1 + 



Am|^ 2 
Am\ 



0.30 (0.44) 



Am^ , 
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With the choices for S, R12 and -R13 made we get using Eqs. ([2]), fl52|) and fl56|) : 



|lm (e'^'^'f/^af^rs)! = C23 C13 (S23C12 + C23S12S 



13 J 



. "32 

sm 

2 



- (C12 + S12S13J 



. "32 

sm 

2 



Arrip \ * .032 

-T — 2 1^12! (■523C12 + C23S12S13) e* 2 - K |i?i3| C23 



1 



ml 



1 + (Cl2 + S12S13) (Ci2 + S12S13 - 2 K COS -^J 



m2 = rrie + 



where k = e^^^^s = sgn(i?i2-Ri3) = ±1 and we have set S23 
|i?i3| specified above and neglected the terms oc s^g in Eqs. 



(82) 
(83) 
(84) 

(85) 



C23 = l/v^5 used I-R12I and 
and 



It follows from Eqs. fl49l) . fl82l) and flHojl that in the case under discussion the following 
relations hold: (032)1 = |e^(27r — 032)!, and frfr {0(32, P23 = 0) = frTri^n — a32, P23 = tt). As 
a consequence we have from Eq. ( 144|) : \YB{a32, P23 = 0)| = \Yb{27i — 032, /323 = 7r)|. Thus, 
we will consider values of 0:32 in the interval [0, 2n] and limit our discussion to (323 = (i.e., 
K = sgn(i?i2-Ri3) = 1). After fixing the value of (323 the only free parameter left in the 
problem being studied is the Majorana phase 032. 

For 032 = 0; 27r, obviously \er\ = and therefore \Yb\ = 0. We can have \Yb\ = also 
even if je,-! 7^ 0, provided the efficiency factor in the expression for \Yb\, Eq. fl44l) . is zero, 
i.e., if the condition given in Eq. (16T1) is fulfilled. For the values of I-R12I and I-R13I considered 
and ^32 having a value in the interval [0, 27r], this condition is satisfied for 0^32 = 1.27r. Thus, 
in this case we have le^-l 7^ but \Yb\ = for 0:32 = 1.27r. 

The absolute maximum of the baryon asymmetry lY^I as a function of 032 is reached 
close to ^32 = vr/2. We have at the absolute maximum for S13 = (0.2): 



1.2 X 10" 



I 0.05 eV 



Ml 



109 GeV 



\Yn 



2.0(2.2) X 10' 




^6) 
^7) 



As we have already indicated, the values of I-R12I and I-R13I we have chosen in this analysis 
maximise not only e^, but also \Yb\. Thus, the observed baryon asymmetry having a value 
in the interval 8.0 x 10"" < \Yb\ ^ 9.2 x 10"", can be reproduced for Mi > 3.6 x 10^° 
GeV. A second (local) maximum of \Yb\ occurs close to 032 = 37r/2 (at 032 = 1.77r). At 
this maximum the value of \Yb\ for S13 = (0.2) is approximately by a factor of 3 (1.8) 
smaller than that at the absolute maximum in Eq. (IHTl) . Given the fact that the fiavour 
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effects of interest can be substantial for values of Mi up to Mi ~ 10^^ GeV, leptogenesis 
can be successful even for rather small values of the Majorana phase 0^2, or more precisely, 
of I sin(a32/2)|. If, for instance. Mi ~ 5 x 10^^ GeV, the observed baryon asymmetry can 
be generated for | sin(a32/2)| = 0.15. 

The results obtained in this subsection are illustrated in Figs. [U and [2] we have already 
discussed, and in Fig. [3] where the dependence of the baryon asymmetry \Yb\ on the 
Majorana phase 032 is shown for Mi ~ 5 x lO^'^ GeV in the case analised above: 6 = 0, 
real R12 and -R13 which maximise \Yb\, i.e., |i?i2| = 0.92 and |i?i3| = 0.39, and P23 = 0, 
i.e., K = sgn(i?i2-Ri3) = 1 (Fig. E]), and P23 = tt, i.e., k = —1 (Fig. H]). The figures have 
been obtained using the best fit values of the oscillation parameters Am^, AtUq, sin^ ^23 
and sin^6'i2. Results for two values of sin ^13 are presented: sin 6*13 = 0; 0.20. Figured 
shows, in particular, that in the case being studied, the predicted baryon asymmetry \Yb\ 
exhibits very weak dependence on sin 6*13 for 032 ^ vr. If sin 6*13 has a value close to the 
existing upper limit, the effect of sin ^13 ~ 0.2 can be noticeable in the region of the local 
maximum of IYbI at 032 = l.Tvr: it can lead to an increase of IYbI by a factor of 1.7. For 
Ml = 10^^ GeV, for instance, we can have successful leptogenesis for 7r/4 ^ 032 ^ 0.97r, 
and, if S13 ~ 0.2, also for lAn ^ 0:32 ^ 1.97r. We note that, as Figs. [H and [2] show, that 
the predicted value of \Yb\ exhibits a relatively strong dependence on the elements of the 
matrix R. 

It follows from the results obtained in the present subsection that as long as | sin(a32/2)| 
is not exceedingly small and Mi ^ 3.5 x 10^° GeV, we can have successful leptogenesis 
even if |si3sin(5| = (Jcp = 0) and the only CP- violating parameter is the low energy 
Majorana phase 032. 



5.2.2 Dirac CP-Violation in ?7pmns and Leptogenesis 

The next question we would like to address is under what conditions we could have a 
successful leptogenesis if the only CP-violating parameter is the Dirac phase in U, i.e., if, 
e^., (323 = 0, vr, the Majorana phase 032 takes a CP-conserving value and sin(a;32/2) = 
Pi . so that 032 = 2Tik, k = 0,1,.... We choose for concreteness again real R12 and R13 
which for 032 = (27r) and /323 = tt (0) maximise \er\ and I^bI, i.e., I-R12P = (1 + 
^yAml/^Aml)-^ = 0.85, |i?i3p = 1 - |i?i2p = 0.15. We will present resuhs for the 
baryon asymmetry \Yb\ also for the values of R12 and -R13 which maximise \Yb\ in the case 
of ^32 = (27r) and /?23 = (tt), |i?i2p = 0.75 and \Rn\^ = 0.25. For 5 = vrg, g = 0, 1, 2, 
the CP-symmetry is not violated and |e,-| = 0. 

For the chosen CP-conserving values of 023 and P23 we get from Eqs. ([2]), ( l52l) and ( l56i) : 

|lm {e'^''U:^Ur3) I = 4 si2 si3 \smS\ , (88) 

^°Obviously, should have a value sufRciently different from the "special" one for which |77(0.6677i^) — 
?/(0.71r^2)| 0. 

By imposing the condition sin(Q!32/2) = wc exclude the possibility of 0132 taking the CP-conserving 
values a32 = tt, Stt, in which case the term associated with the violation of CP-symmetry due to the matrix 
R will dominate in |Im([/*2C/r3)|, and correspondingly in \er\- If P23 ~ 7''/2,37r/2, the indicated possibility 
would be avoided if a32 takes a CP-conserving value and | sin(Q!32/2)| = 1. 
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rrir = a/ Am^ 



A 



-^^J I-R12I (S23C12 + C23S12S13 e'^) - k' |i?i3| C23 Ci3 



""^^ . |(c23Cl3-K-S23Cl2)+C233l2Sl3e-^r k' = e^(^- + "i") = ±1 ,(89) 



1 



i^Am^ 



where in Eq. ( 1891) we have used the values of I-R12I and I-R13I specified earlier. Taking into 
account that |Im(e*^23f/*^f/^3)| ^ 0.027|si3 sin we find from Eq. ([79]): 



' I - 2 2 X 10-^ [ VM ] \ 



(90) 



Thus, the maximal asymmetry le^-l in the case of CP- violation due only to the Dirac phase 
S in ?7pMNS is approximately by a factor of 6 smaller than the maximal asymmetry due to 
violation of the CP-symmetry by the Majorana phase 0^32 of ?7pmns- 

Given r7v, one can determine from Eq. (1851) . It follows from Eqs. (1891) and (I90p 
that \Yb{S)\ = \Yb{2tt - 5)\. One has to analyze the cases of k' = e*(°32/2+fe) ^ 
(032/2 + P23 = 2TTk, k = 0,1,...) and k' = -1 (032/2 + P23 = vr(2A; + 1), A; = 0, 1, 2, ...) 
separately. 

For k' = —1 we find irfr = 0.25y^AmJ > m2 = O.OSy^AmJ. For the values of the 
parameters employed in this analysis (chosen, in particular, to maximise |e^| and 11^1); ^ 
and m2 exhibit weak dependence on sin ^13 (and therefore on 6), which can be neglected. 
This implies that the maximum of the baryon asymmetry 11^1 as a function of the Dirac 
phase S, will take place at values of 5 = (7r/2)(2fc + 1), k = 0,1,..., for which \er\ also 
has a maximum. Using Eq. ( l45i) to calculate the relevant efficiency factors r]{0.66rrfr) and 
77(0.717712), we get from Eq. (jH]): 

^ o s / ]\/r \ 

(91) 



|y.|.2.8xl0-|.n, (a)(To^) 



The asymmetry of interest is predominantly in the lepton number Lg + L^. Thus, in order 
to reproduce the observed baryon asymmetry, taken to lie in the interval 8.0 x 10~^^ ^ 
\Yb\ ^ 9.2 X 10~^^, S13I sm6\ and Mi, in the case anahsed, should satisfy 

Given that Si3|sin5| ^ 0.2, the lower bound in this inequality can be satisfied only for 
Ml ^ 2.9 X 10^^ GeV. Recalling that the fiavour effects in leptogenesis of interest are fully 
developed for Mi ^ 5 x 10^^ GeV, we obtain a lower bound on the values of |si3sin5| and 
Si3 for which we can have successful leptogenesis in the case considered: 

Isin^ia sin5| > 0.11, sinei3>0.11. (93) 
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The lower limit fl93!) corresponds to 

|Jcp| ^ 2.4 X 10-^ (94) 

where we have used the best fit values of sin 2^12 and sin 2^23- Values of S13 in the range 
given in Eq. fl93p can be probed in the forthcoming Double CHOOZ [HH] and future reactor 
neutrino experiments . CP- violation effects with magnitude determined by | Jqp \ satisfy- 
ing are within the sensitivity of the next generation of neutrino oscillation experiments, 
designed to search for CP- or T- symmetry violations in the oscillations [67] . Actually, since 
in the case under discussion the wash-out factor |?7b| = |?7(0.66rrv) — ?7(0.71m2)| in the ex- 
pression for \Yb\ practically does not depend on S13 and 6, while both \Yb\ oc |si3sin5| 
and |Jcp| oc |si3sin(5|, there is a direct relation between \Yb\ and |Jcp| for given neutrino 
oscillation parameters, -R12, -R13 and Mi: 



\Y, 



B\ 



Mi/(10ii GeV) 



3.0 X 10"*^ \r]B\ I JcpI = 1-3 x 10"^ | Jcp| , (95) 



where we have used the best fit values of the neutrino oscillation parameters, I-R12I = 0.92, 
|i?i3| = 0.39 and k' = -1. 

For n' = +1, the maximum of \Yb\ as a function of I-R12I (I-Risl) takes place in the 
case being investigated at I-R12I = 0.86 (I-R13I = 0.51) and we will employ this value of 
I-R12I (I-Risl) in the remaining part of the current analysis. If = -|-1, a strong cancellation 
between the terms in the round bracket in Eq. ( l89i) takes place and we have typically irfr -C 
rn2- Indeed, for S23 = C23 = l/"\/2, S12 = 0.2 and 6 = n/2 one finds irfr — 2.5 x 10'"^ \/ Am\ 
and m2 = 0.28^/ Am\. Using Eq. ( l45i) to calculate the corresponding efficiency factors, we 
get from Eq. (1471) for the baryon asymmetry: 

IFbI =3.6 (0.38) X 10-i3|gi^^| (^-^^^^ . (96) 

The numerical factor in the round brackets corresponds to the value of I-R12I = 0.92 (I-R13I = 
0.39) which maximises |e^|, but does not maximise \Yb\. 

The result we get for the maximal value of \Yb\, i.e., for |_Ri2| = 0.86 (I-R13I = 0.51), are 
quite similar to those we have obtained in the case of k' = —1, the only difference being 
that now the generated lepton asymmetry is predominantly in the tau lepton charge. The 
observed baryon asymmetry lying in the interval 8.0 x 10~^^ ^ \Yb\ ^ 9.2 x 10~^^, can 
be reproduced if Mi ^ 2.2 x 10^^ GeV. For the flavour effects to fully develop one must 
have Ml ^ 5 x 10^^ GeV, which, together with the requirement of successful leptogenesis, 
imphes 

I sin ^13 smS\> 0.09 , sin ^13 > 0.09 . (97) 
This lower limit corresponds to 

|Jcp| > 2.0 X 10~^ (98) 

where we have used again the best flt values of sin 2^12 and sin 2^23- The ranges of values 
of sin ^13 and of |Jcp| we flnd in the case being studied are also within the sensitivity 
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respectively of the planned ^13 reactor neutrino experiments [651 ES] and of the neutrino 
oscillation experiments on CP- (T-) violation [U7] . 

It is interesting to note that, as it follows from Eq. flUUl) . for |i?i2| = 0.92 and I-R13I = 
0.39, which maximise \er \ but do not maximise \Yb\, the predicted baryon asymmetry \Yb\ 
is smaller by approximately one order of magnitude in spite of the fact that the indicated 
values of |i?i2| and |i?i3| are rather close to the values I-R12I = 0.86, |i?i3| = 0.51 which 
maximise \Yb\. As |i?i2| increases beyond 0.86, 11^1 rapidly decreases, as is also clearly 
seen in Fig. [TJ Actually, for |i?i2| = 0.92 and I-R13I = 0.39, the observed baryon asymmetry 
cannot be reproduced for Mi < 10^^ GeV and sin ^13 < 0.20: we get \Yb\ ^ 4 x 10~^^ 
Figures [U and [2] show also that for real R12 and -R13, satisfying I-R12P+ l-RisP = 1, and CP- 
violation generated only by the Dirac phase in t/pMNS, we can have successful leptogenesis if 
both sin ^13 and I-R12I have relatively large values: the observed baryon asymmetry cannot 
be reproduced if I-R12I ^ 0.6 and/or sin6'i3 ^ 0.1. 

The results obtained in the present subsection are illustrated in Figs. [T] and [21 as well 
in Fig. In Fig. [3] we show the dependence of \Yb\ on the Dirac phase 6 which was varied 
in the interval [0, 27r], for S13 = 0.1; 0.2, Mi = 5 x 10^^ GeV and for k' = +1 and k' = -1. 
The correlation between the rephasing invariant Jcp which controls the magnitude of the 
CP-violation effects in neutrino oscillations and the baryon asymmetry Yb is illustrated in 
Fig. [HI for Si3 = 0.2, Mi = 5 x 10^^ GeV and k' = +1. Both figures are obtained for real 
R12 and i?i3 which maximise \Yb\. 

In conclusion of the present subsection we note that one can treat in a similar way the 
alternative possibility of P23 = 7r/2 (37r/2). In this case we have oc Re(U*2Ur3)- The 
CP symmetry will be violated at low energies if both Re(f/*2f/r3) 7^ and lm(U*2Ur3) 7^ 0. 
Maximal asymmetry le^l is obtained for cos(a;32/2) = ±1, and, if S13 is non-negligible, for 
cos((5 — 032/2) = ±1. These conditions are satisfied for the CP-conserving values of 032 and 
6: a-s2 = 0; 2tt, 6 = 0; n. The CP-symmetry is broken by the matrix R {(323 = 7r/2 (37r/2)). 
It is easy to convince oneself that the expression for the asymmetry |e;|, I = e,fi,T, for 
l323 = 7r/2 or 37r/2, 032 = (27r) and S = 0{ n), coincides with the expression for the 
same asymmetry in the case respectively of /323 = or tt, 032 = vr (37r/2) and 6 = {2n). 
Although, for /323 = vr/2; 37r/2 we also have max(r) = 1.2, the maximum of r corresponds 
to |i?i2p = = 1.22, {Ri^l^ = {Rul^-l = 0.22. Therefore the wash- 

out factors mi will differ from those in the case of /523 = 0; vr and maximal asymmetry |e^|. 
The values of I-R12I and I-R13I which maximise le,-! are not guaranteed to maximise also the 
baryon asymmetry \Yb\. Further investigation of this case is, however, outside the scope of 
the present study. 

5.3 Inverted Hierarchical Light Neutrino Mass Spectrum 

Given the inequalities m-^ <^ mi < m2, it seems natural to assume that the terms oc ^/rf^ 
are sufficiently small and give negligible contributions to the lepton flavour asymmetries e/ 
and to the wash-out mass parameters mi. In the "one- flavour" approximation this case was 
studied in [28]. It was found that the lepton asymmetry is strongly suppressed by the factor 
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AmQ/Am^ and leptogenesis can produce the observed value of the baryon asymmetry only 
if the lightest heavy Majorana neutrino Ni has a relatively large mass: Mi ^ 7 x 10^^ GeV. 
Using Eq. fjlUj) and the fact that in the case of interest we have mi 2 = a/ Am^ and 
(m2 — mi) = AttIq/ (2-\/Am^), it is not difficult to obtain the conditions under which 
the indicated terms oc y/m^ in ei can be neglected. The latter depend on whether the 

R11R12 = e*^^^|i?ii-Ri2| is real or purely imaginary. In the case of P12 = Pn + Pu = vrg, 
g = 0, 1,2, we get: 



i „ 3 

mg \ /Am^V |i?i3| 



^Am^ / KAmlJ |/2i2(ii 



<1, /5i2 = vrg, g = 0,1,2,.... (99) 



Since 2(Am^/Am|)3 = 26.4 ^ 1, this inequality suggests two possibilities. 

i) Equation ( 199|) holds and the terms oc y^ms in ei and m/ are indeed negligible. The 
simplest realisation of this possibility corresponds to [28] setting R13 = 0, which in turn 
implies the decoupling of the heaviest RH Majorana neutrino N^. 

ii) The alternative possibility is that terms oc ^/rr^ in e/ and rui are dominant in spite of 
the fact that ma <^ mi, m2. This would require the ratio in the left-hand side in Eq. (|99l) 
to be much bigger than 1. A possible simple realisation corresponds to setting i?ii = or 
i?i2 = 0. Since in the latter case le^l oc ^/ms/rri2, the asymmetry \ei\ will not be suppressed 
only if m3 is sufficiently large. The latter condition is satisfied for values of in the 
interval 10~^ y^Amg ^ ^ O.b^y ArriQ, for which we still have m^, <C mi 2. 

If, however, P12 = n/2{2q + 1), q = 0, 1, 2, i.e., if RuRi^ = ±i I-R11-R12I, we obtain a 
very different condition: 

< 1 , Pu = (vr/2)(2g + 1), g = 0, 1, 2, ... . (100) 





2 










-Rl2(ll) 



For m3 -C ^7^12,, this condition can be naturally satisfied if I-Risl is sufficiently small and, 
in particular, if |i?i3| = 0. 



5.3.1 The case of Real RuRu and A^3 Decoupling {R13 = 0) 
The terms oc ^/rf^ in e; and m/ are negligible and we get: 



/ Am|\ / Am|\ ^ |i?ni^i2| . , ^ _ 



(101) 

where j3i2 = aiglRuRu) ■ The phase /3i2 parametrises the effect of CP- violation due to the 
matrix R in the asymmetries e^. In the case of CP-invariance, /5i2 = or 7r/2 depending on 
whether Im([/f^[//2) = or Re(f/fif/;2) = 0, respectively, and e; = 0, / = e, /i, r. Note that 
the asymmetries ei are suppressed by the factor AmQ/(2Am^) = 1.6 x 10~^ with respect 
to the asymmetries e/ we have obtained in the case of NH light neutrino mass spectrum. 
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From Eq. ([2]) we get for the quantity lm{e^^^^U*iUr2) of interest: 

Im {e'^^'U:,Ur2) =-Im [e'^^' e'^ (012^23 + S12C23S13 e^') (512^23 - C12C23S13 e"'') 

"(102) 

where 021 is the Majorana phase which enters also into the expression for in (/5/3)ov- 
decay. The wash-out mass parameters mi for the three lepton asymmetries read [28] : 



mi^^Aml\RuU;^ + Ri2U;2\, l = e,iJ,T. (103) 
Using the unitarity of the PMNS matrix U we obtain: 



me + m^ + m, = ^ Ami (1^11^ + 1^12^) • (104) 

In what follows we limit our discussion to the case of real -Rn and R12. In this case we 
have I-Riip + I-R12P = 1, and Pu = 0; vr which correspond to sgn(i?iii?i2) = +1; (— 1)- 
Obviously, the values of |-Rii| = I-R12I = l/\/2 maximise the asymmetries e/. As can be 
shown, for generic values of the phases a2i and 6, these values of and |i?i2| maximise 
also \Yb\- Now we have 



me + m^ + mr = y Am\ . (105) 

Using the best values of the neutrino oscillation parameters sin^ 26'23, sin^ 6*12, Am^ and 
Am\, given in Eq. ([5]) and the upper limit S13 < 0.2, we get maximal le^l and IYbI as can 
be shown, for 0^21 = vr/2, 6 = n and P12 = 0: 



Clearly, for Mi ^ 10^^ GeV for which the flavour effects in leptogenesis can be substantial, 
it is impossible to reproduce the observed baryon asymmetry in the case of IH light neutrino 
mass spectrum, R13 = and real Rn and R12, l-RiiP + I-R12P = 1- The main reason for 
this result lies in the fact that the lepton asymmetries ei are suppressed by the factor 
Am|/(2Am^) = 1.6 x 10~^, while the wash-out effects are rather large owning to the 
constraint f llOSp . 



5.3.2 Generating Yb Compatible with the Observations 

Consider next the case of = (7r/2)(2g+ 1), q = 0,1,2,.... Now the product R11R12 
is purely imaginary: R11R12 = m|-Rii-Ri2|, in = e^^^^, k = ±1. We shall assume further 
that condition fllOOl) holds and thus the terms oc y/m^ in ei and mi are negligible. A simple 
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realisation of this scenario corresponds to -R13 = {N^ decoupling), and to — |i?i2| 

1. Under the indicated conditions we have: 





2 _Rii_Ri2 






Ru 


2 + 


-R12 


2 



|Re {U*^U-r2)\ , I = e,iJ,T, 



:io8) 



where 

|Re {U:,Ur2)\ 



C12S12-S23 COS ■ 



"21 



C23S23S13 



COS 



/a2i _ 
V 2 



cos 



(109) 

and we have neglected the term oc 5^3012512023 in the expression for |Re(f/*^?7r2)|- The 



maximal value of the factor r' = 2\RiiR 



12 



R 



11 



\Ri2?) is r' 



1 and is reached for 



l-RiiP ^ 1- However, even for I-Rup = 1.5 and I-R12P = 0.5, r' is rather close to its 
maximal value: r' = V3/2 = 0.87. In the case of purely Majorana or Dirac CP-violation 
from the PMNS matrix we obtain: 



\Re {U*^Ur2)\ = S23 (C12S12S23 ±Si3C23 cos 26^12) 

\Re {U*^Ur2)\ = S13C23S23 |sin5| , cos ^ = . 



cos ■ 
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6 = 0,71,2-11, (110) 
(111) 



A. Majorana CP- Violation from U] 



PMNS 



In what follows we first set S13 = 0. The CP-symmetry is violated by the Majorana 
phase ^21 only. Both \er\ and \Yb\ vanish for 021 = 7r(2g + 1), q = 0, 1, .... The wash-out 
parameters r?v and are given by 



= \/Am\ 



2 2|Dl2,2 2|Dl2,r) 2|DDl - '^21 

'12 ■^23 l-fLiil + C12 523 |-n:i2| + 2KC12S12S23 |-n:ii-n:i2| sm — 



:ii2) 



m2 = ^Aml{\Rii\^+\Ru\^) -m,. (113) 

Obviously, we have |eT-(a2i)| = |e^(27r — 021) |, rrv(a;2i) = mT-{2n — 021), and therefore 
|^b(q^2i)| = |^B(2'/r — 021) I . Numerical calculations we have performed show that the 
maximal value of the baryon asymmetry \Yb\ in the case under discussion corresponds to 
l-Riip = 1.1 — 1.4 (see Fig. [7j) and we will use these values of I-Rup and the corresponding 
values of I-R12P = 0.1 — 0.4, in our further analysis. The results we obtain for \Yb\ depend 
strongly on whether K = +lorK = — 1. 

For K = +1 and I-Rup = 1.1 — 1.4, one has rrv,?^ ^ 2 x 10~^ eV, which corresponds 
to a strong wash-out regime. The baryon asymmetry is maximal for a2i = 27rg, g = 0, 1, 
for which both |e^| and the efficiency factor |?7b| = |?7(0.66r7v) — ?7(0.71m2)| are maximal PI 
(I cos(a2i/2)| = 1). Numerical studies show that the absolute maximum of IYbI for k = +1 
is reached at |i?iiP = 1.1 (Figs. [7] and El). 



-^^Note that for a2i = 27r(7, q = 0,1, and S13 = we have l-m{U*iUr2) = and Jcp = 0. Corre- 
spondingly, the CP-invariance is not violated by a2i and 5; it is violated at high energies by the matrix 
R. 
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In the case of k = — 1 and values of |-RiiP indicated above, which maximise 1^51; there 
is a strong partial compensation between the three terms in m^. As a consequence, the 
efficiency factor \riB\ and the asymmetry \Yb\ are for, e.g., 021 = vr/2 and the values of 
which maximise \Yb\, approximately by a factor of 5 bigger than the values they 
have in the case of k = +1. Moreover, the maximum of \Yb\ takes place at I-Rup = 1.4 and 
021 — 27r/3; 47r/3 (Figs. [7] and [9]), rather than at I-Rup = 1.1 and 021 = 27rg, g = 0, 1, .... 
In both cases of k = +1 and k = — 1 we have \Yb\ = for ^21 = 7r(2g + 1). 

We get for the maximal baryon asymmetry (i.e., at I-Rup = 1.1 and 0:21 = 0; 2tt for 
K = +1, and at = 1.4 and 021 = 27r/3; 47r/3 for k = -1, see Figs. Ul-E)- 

WM - 1.5 (0.5) X 10-'^ (^) (J^) . . ^ -1 (. ^ . (114) 

It follows from the preceding analysis that the observed value of the baryon asymmetry 
IIbI ^ 8 X 10~^^ can be reproduced in the case being studied for Mi ^ 5.3 x 10^*^ GeV if 
K = —1, and for Mi ^ 1.6 x 10^^ GeV when k = +1. Since both the baryon asymmetry 
\Yb\ and the effective Majorana mass in (/3/3)oiy-decay, \ {m) \ depend on the Majorana phase 
021, for given values of the other parameters there exists a direct correlation between the 
values of \Yb\ and \{m)\. The latter is illustrated in Fig. [TUl 

B. Dirac CP- Violation from t/pMNS (<^2i = t^) 

One can treat in a similar manner the case of S13 7^ and CP-violation generated only 
by the Dirac phase 5 in f/pMNS, ^21 = 7r(2A; + 1), k = 0, 1, .... We have e*"^^/^ = in', where 
= 1 for A; = 0, 2, 4, and = -1 if A; = 1, 3, .... It follows from eqs. ffTOSD and ffTTOj) 
that in this case the asymmetry \er\ oc |si3 sm6\ and one can expect the baryon asymmetry 
\Yb\ to be suppressed by S13. Obviously, we have \er\ = and IYbI = for 5 = 0; vr; 27i. 
For the wash-out mass parameters r?v and rn2 we get: 

m; ^ \l |S23 (Si2 \Rll\ - Kk'ci2 \Ri2\) - C23Si3e~'^ (Ci2 \Rll\ + Kk' S12 |i?12|)f , 

_ _ (115) 

and m2 = \/ Am^(|i?iip + I-R12P) — rrir. The two possibilities kk' = —1 and kk' = +1 lead 
to drastically different results. 

For kk' = — 1 , the term with the factor S23 in the expression for in Eq. (11151) gives the 
dominant contribution and determines the magnitude of r?v- For, e.g., 5 = tx/2 (for which 
|e.r| is maximal), the maximum of the baryon asymmetry \Yb\ takes place at |-Rii| = 1.05 
(I-R12I — 0.32). For these values of 5 and |-Rii|, both wash-out mass parameters satisfy 
rn^2 ~ 2 X 10~^ eV. Thus, the baryon asymmetry is generated in the strong wash-out 
regime. Correspondingly, the efficiency factor |?7b| is relatively small, |?7b| ^ 6 x 10~^. The 
observed value of the baryon asymmetry \Yb\ = (8.0 — 9.2) x lO"^'^ can be reproduced 
only if si3 = 0.2 and Mi = (4 - 5) x 10" GeV; for S13 = 0.1, this requires a value of 
Ml = 8 X 10" GeV. 

We obtain completely different results in the case of kk' = +1. In this case there can 
be a deep mutual compensation between the two terms in the bracket multiplied by S23 
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in the right-hand side of Eq. (11151) and we can have rrv — (1-5 — 2.0) x 10~^ eV for the 
values of |-Rii| (I-R12I) which maximise the baryon asymmetry \Yb\. Correspondingly, \Yb\ 
can be generated in the weak wash-out regime and can be much larger than in the case of 
kk' = —1. 

More specifically, the maximum of lYsl with respect to the Dirac phase 6 and |-Rii| 
(|i?i2|) takes place approximately at (or relatively close to) 6 = 7i/2; 37r/2 for any S13 ^ 0.2 
of interest, and at |i?ii| ^ 1.30; 1.60 (|i?i2| = 0.83; 1.25) for S13 = 0.20; 0.10, respectively 
(Figs. [II]and[I2D; for S13 < 0.02 it is located at \Rn\ ^ 1.07 (|i?i2| = 0.38). We have 
used again the best fit values of Am\, sin 2^12 and sin 2^23 to obtain the location of the 
maxima in |-Rii|. In the case of S13 = 0.10 there exists a second local maximum of \Yb\ at 
|i?ii| = 1.15 (I-R12I — 0.57), at which \Yb\ is approximately by a factor of 1.4 smaller than 
at the absolute maximum at |-Rii| = 1.60 (see Fig. [TT]). The positions of the indicated 
maxima of \Yb\ are determined essentially by the position of the absolute maximum of the 
efficiency factor |?7b| = |?7(0.66r7v) — ?7(0.71fra^)|. The latter corresponds approximately to 
rrv = 1.5 — 1.8 x 10"^ eV and negligiblecl |?7(0.71m^)|. At the maximum, Ir^^l = 6.7x 10~^. 
In the case of S13 = 0.20, as can be shown, irfr ^ 2.3 x 10^"^ eV and the minimal value of 
frfr corresponds to negligible S23('Si2|-Rii| — Ci2|-Ri2|)^ — 0. The latter condition is fulfilled 
approximately for \Rii\^ = cl^/ cos 29 12 = 1.75 and |i?i2p = 5^2/ cos 26^12 = 0.75, which 
is very close to the value I-Rup = 1.69 obtained by numerical calculations. Thus, for 
Si3 = 0.20 we have for which maximises |?7b ^23-513(^12 1 -Rill + 

Si2\Ri2y = (c^gs^g/ COS 2^12) VAmi. 

In the case of Si3 = 0.10, both terms in Eq. (I115p "conspire" to produce r7v = 1.75x10^^ 
eV, and thus maximal \r]B\ (and llsl), at \Rii = 1.60. The second local maximum of \Yb\ at 
l-Riil = 1.15 corresponds to irfr — 10~^ eV. Again both terms in Eq. f lll5p contribute, the 
term oc c^^sl^ being approximately by a factor 1.5 smaller than the term oc S23. In contrast, 
the local minimum of \Yb\ at |-Rii| = 1.28 (Fig. [TT]) is associated with ttv = 6 x 10~^ eV. 
In this case the contribution of the term oc S23 in is negligible. 

Finally, for S13 ^ 0.02, the term with the factor oc C23S13 in the expression for rTv, 
Eq. f lll5p . plays no role in the determination of the maxima of \t]b\ and \Yb\. Thus, in this 
case, in particular, \Yb\ depends on S13 and 6 only through and we have \Yb\ oc S13I smS\. 

For the maximal value of the baryon asymmetry \Yb\ for S13 = 0.20 (0.10) and 6 = 7r/2 
(|i?ii| ^ 1.30 (1.60)) we get (Fig. M>: 

1^,1 . ,r ,1.0) X 10-^ (^) (^) . * ^ ^ , .3 ^ 0.20 ,0.10) . ,116) 

It follows from the above result that if, e.g., S13 = 0.10, we can obtain baryon asymmetry 
\Yb\ ^ 8 X 10~^^ compatible with the observations for Mi ^ 8 x 10^° GeV. 

What is the minimal value of S13I sin5| for which we can have successful leptogenesis? 
For sufficiently small values of S13 we get maximal |?7b| = 6.5 x 10~^ {frfr = 1.75 x 10~^ eV) 
at = 1.07 independently of the value of S13 (see the discussion preceding Eq. (11160 ). 

^^Indeed, since + |i?i2p > 1, wc have > 5.0 x 10^^ eV and therefore |?7(0.71to^)| < 2 x 10"^. 
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The baryon asymmetry at the maximum is given by: 



,,,.3.x:o-.^.3|s...|(S)(^). (n. 

Thus, for Ml ^ 5 x 10^^ GeV, we can have successful leptogenesis and get \Yb\ = (8.0 — 
9.2) X 10"^^ provided 

Isin^ia sin(5| > 0.02, sin^i3>0.02. (118) 

The preceding lower limit corresponds to 

|Jcp| ^4.6 X 10-^ (119) 

Values of sin6'i3 and of |Jcp| as small as 0.02 and 4.6 x 10~^, respectively, can be probed 
in neutrino oscillation experiments at neutrino factories [67] . 

The results in the case under discussion are illustrated in Figs. [H] and [121 In Fig. 
[TT] we show the baryon asymmetry as a function of |-Rii| for 6 = 7r/2, while Fig. [T2] 
exhibits the dependence of \Yb\ on the Dirac phase 6 for the values of |-Rii| from the in- 
terval l-Riil = (1.05 — 1.7), which maximise \Yb\- The results presented in both figures 
are for sis = 0.1; 0.2, 021 = tt {k' = +1), k = +1, k = -1 and Mi = 2 x 10^^ GeV. 
In Fig. [13] we show the correlation between the rephasing invariant Jcp (in blue) and the 
baryon asymmetry Yb in the case under discussion, (CP-violation due to the Dirac phase 
S = tt) for = 0.2, Ml = 2 X 10^^ GeV, k = +1 and |i?ii| = 1-3. The Dirac phase 6 
is varied in the interval [0,27r]. 

C. CP-Violation due the Dirac Phase in t/pMNS ^i^d R (a2i = 0; 27rj 

We will consider next briefly the "mixed" case of CP-violation corresponding to a2i = 
27rk, k = 0,1,..., S taking values in the interval [0,27r], and purely imaginary R11R12, 
-R11-R12 = "ii^lRiiRul, ft; ± 1. For 6 = 0; tt, this case provides an interesting example of 
breaking of CP-symmetry due to R11R12 being purely imaginary. 

It proves convenient to write e*"^^/^ = k', where k' = ±1. We have in the case under 
discussion: 

|Re {U*iUr2)\ = \ci2 S12 S23 - C23 S23 cos 26*12 si3 cos5| . (120) 

Obviously, the term cx Si3Cos5 plays a sub-dominant role in the asymmetry \eT-\. The 
wash-out mass parameter ??v reads: 

mr = y^Ami [5^3 (s?2 \Rii\^ + 4 \Ri2\^) + C23 sis (4 l^nP + ^2 1^12^ 

- 2C23S23S13 (C12S12 cos5 + kk' |i?iii?i2| sin(5)] . (121) 

It follows from eqs. ffTUH]) . ([nUD and ([HID that |e^(5)| = |e^(27r-(5)|, r^(5, kk') =r^(27r- 
6, —k,k'), and therefore \Yb{S, kk,') \ = Yb{1t{ — 5, —kk')\. 
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For ^ 0.01, the terms oc S13 and oc sl^ in the expressions for |e,-| and r?v can be 
neglected. The CP- violation effects are practically due to the purely imaginary R11R12: 
the baryon asymmetry Yg does not depend on S13 and S. However, even in this case it is 
possible to have successful leptogenesis. The baryon asymmetry Yb can be generated in a 
regime which is intermediate between the weak and strong wash-out ones: irfr ^ 8 x 10^"^ 
eV. The asymmetry \Yb\ has two very similar (in magnitude) maxima at |-Rii| = 1.05; 1.12. 
At these maxima we obtain: 

|y.| ^ 5^2 (5,1) X 10- (j^) [Jl^) . ..3 = 0. \Rn\ ^ 1.05 (U2) . (122) 

Thus, one can have \Yb\ = (8.0 - 9.2) x 10"" for Mi > 1.6 x 10^^ GeV. 

If Si3 ^ 0.01, the effects of the Dirac CP-violating phase 6 are non-negligible; for 
Si3 = (0.1 — 0.2), they are substantial for any value of 6. The magnitude of the baryon 
asymmetry \Yb\ depends on the value of S13 and, through the wash-out mass parameters 
rrv and m2, on the sign factor kk,'. In the case of 5 = 7r/2, the maximum of \Yb\ is located 
approximately at |-Rii| = (1.05 — 1.10) (Fig. [T^ . For kk' = +1, the baryon asymmetry at 
the maximum is by a factor ~ 2 larger than that for kk' = —1. We can have successful 
leptogenesis in the case of S13 = 0.1 (0.2) for Mi ^ 1.2 (1.0) x 10^^ GeV. These results are 
illustrated in Figs. [TH and [I5l In Fig. [16] we show the correlation between the rephasing 
invariant Jcp and the baryon asymmetry 1b for 021 = 0, S13 = 0.2, n = +1, \Rn\ = 1-05 
and Ml = 2 x 10^^ GeV. The Dirac phase 6 is varied in the interval [0,27r]. 



As we have indicated at the beginning of this subsection, the observed value of the 
baryon asymmetry can be reproduced in the case of real Ry, j = 1,2,3, if the terms 
cx ^/m^ in e; and are dominant in spite of the fact that rris <^ mi,m2. A simple 
realisation of this possibility corresponds to having Ru = or R12 = 0, and sufficiently 
large 777,3 still obeying the inequality <^ mi^2- The latter conditions can be satisfied for 
m3 possessing a value in the interval 10~'^ AtUq ^ 7713 ^ 0.5^/Am^. A more detailed 
discussion of this case will be presented elsewhere. 



5.4 Quasi-Degenerate Light Neutrino Mass Spectrum 



We turn now to the QD spectrum, for which mi ~ ?7i2 — fn^ ^ y^Am^, \/ Am^ and we 
take as conventional lower limit on the masses m\ ^0.1 eV. In this case the contribution 
of all the masses to the CP-asymmetry needs to be taken into account. We study 
the two cases in which the products of the parameters in R entering in the asymrnetry, 
RiiRij, are either real or purely imaginary. This corresponds to having Pij = Pu + [3ij = 
2/c7r/2, {2k + l)7r/2, k = 0,1, respectively. We consider first the case of RuRij real and 
equal to ±\RuRij\. The CP-asymmetry can be written as: 



3Mi 



167rw2 \R 



li\ 



1712 — TTT'l 



\RiiRi2\e"^'nm(u;^Ur, 



+ {ms-m,)\RuRi3\e^'''nm(u:,Ur3) + ims-m2)\Ri2Ri3\e'^^^^^ , (123) 
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3Mi 1 



Am, 
2m 



2 . s 



2mi 



(124) 



where we have neglected terms of order Amg/Am^. The signs ± in Eq. f ll24p refer to the 
quasi-degenerate spectrum with a normal or inverted hierarchy. This information could 
not be obtained in experiments which are sensitive to the overall neutrino mass scale, as 
neutrinoless double beta decay or direct neutrino mass searches. However, by exploiting 
matter effects, long baseline neutrino oscillation and atmospheric neutrino experiments 
might be able to establish if the spectrum is with normal or inverted hierarchy. Neglecting 
the terms proportional to Amg, we can further simplify the expression in Eq. (11241) : 



SMimi 1 Am\ 

T I p.. 12 9^2 *23C23Ci3 1 Rl3 \ 



X 



' IGttv^ 2m2 

lu if ■ "31 C23 . /«31 if • "32 C23 . /a32 . 

|i?ii |(si2 sm — -ci2— Si3 sm ± |i?i2 1 (- Ci2 sm — - si2— Si3 sm [-^-o 



(,125) 



where ± refer to the case of Pu = P12 and P12 = kn + Pu, k = 0,1,..., respec- 
tively. Notice that, in general, the asymmetry is suppressed by ^m\/m\. For compar- 
ison with the NH case studied in Section 15.11 we can have a mild suppression of order 
(Am^)'^/'^/(mi(Am|)^/^) ~ 0.2 (1 eV/mi) for large values of mi. If i?i3 is negligible {N^ 
decoupling), the dominant contribution is proportional to Amg which amounts to an ad- 
ditional suppression factor of Amg/Am^ ~ 0.03 [2H]- Typically, we get a CP-asymmetry 
of order 

Ml O.leV, 

~ 1.2x10-^— -i— i?i3i?ii 

10^^ GeV mi 

0.55 sin ^ - 0.17^ sin (^ - 5) - 0.84 sin ^ - 0.11^ sin f ^ - 6] (126) 
2 0.2 ^ 2 ^ 2 0.2 V 2 / ^ ^ 

where we have taken for definiteness i?i2 = -Rn and \Rii\^ = 1. The asymmetry 
decreases linearly with mi and we have evaluated it for the minimal value of mi which is 
allowed for the QD spectrum. Notice that, as far as 031, a32 is not too small, sin ^ ^ 
0.1 (S13/O.2), the contribution of the 6 phase will be subdominant. The phases 021 and 031 
enter in the effective Majorana mass parameter for neutrinoless double beta decay. For the 
QD light neutrino mass spectrum, we have 

1(^)1 = ^1 |c?2C?3 + 44^^"" + . (127) 

The dependence on the phase difference (031 — 25) is suppressed by sfg and can be ne- 
glected. In principle, the phase 021 can be measured in future neutrinoless double beta 
decay experiments if a sufficient precision will be achieved in the measurement of | (m) | and 
of the masses, and if the problem of the computation of nuclear matrix elements will be 
solved [58]. The prospects of measuring the other Majorana phase a^i are far beyond the 
sensitivities of the present and future planned experiments. 
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In the case of RuRij purely imaginary, i.e. Pu+Pij = (2/c+l)7r/2, = 0, 1, in general, 
there is no suppression as the masses enter in the CP-asymmetry via the combination 
nii + nij. We get 



3Mimi 



1 



EAR 



li\ 



X 



± |i?n^l2|Re(f/;if/.2) ± |i?lli?13|Re(t/;if/.3) ± \Rl2Rl3\MU:2Ur3: 



:i28) 



where the +(— ) refer to = 7r/2(37r/2). The expression in Eq. 
but can be simplified if we neglect ^13. In this case it reads: 



is rather lengthy 
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±|i?iii?i2|(-Si2Ci2S^3 cos^) ± |i?iii?i3| (S12S23C23 ^^^-f-) =t 1^12^13 1 (-C12S23C23 COS ^|29) 



As we have neglected the terms oc S13, there is no dependence on the phase 6. However, 
both Majorana phases enter in the expression for e^. 

Let's now turn to the wash-out factors. Using the unitarity condition on U, we find: 



7712 + 



mi \R 

k 



lk\ 



(130) 



In the case of real Rik, for instance, this implies [2H] + ^ = nT'i- Therefore, we can 
expect that the wash-out mass parameters will typically be much larger than 3 x 10~^ eV, 
leading to a strong suppression of the baryon asymmetry. More specifically, r?v is given 
by: 



mi \RnU:, + RuU:^ + RnU:, 



(131) 



We will study first the case of real Rij. Taking, for example, Rn = R12 = R13, we 
obtain rTv = 0.053 mi, for 0^21 = tt and asi = n. In this case the strong wash-out regime 
formulas apply and we have 



7/(390/589 r^) = 3.3 x 10"^ 



0.1 eV\i-i6 



mi 



(132) 



resulting in a baryon asymmetry which is substantially smaller than the observed one. 

A larger efficiency factor can be achieved in the regime interpolating between the ones 
of strong and weak wash-out effects. However, it should be noticed that ?7(m^) and 
depend on the same parameters and they cannot be maximized independently. In fact, for 
real Rij, we can rewrite the CP-asymmetry as: 



3Mimi 



Ami 
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327r^;2 7nlJ2^\Rl^? 
SMinii Ami 



ImiRrsUrsY^RiP^r^ 
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(133) 
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Therefore, while the efficiency factor rj{rrfT-) increases when frfr decreases as the wash-out 
regime changes from strong to weak, le,-! decreases accordingly. In the strong wash-out 
regime in the r flavour and r?v <^ '^i, the baryon asymmetry has an upper bound: 



-3 conX 



12 3MimiAmi,„ , /m; / 0.2 x 10"^ eV 589 



- 37^ 32vr^;2 1^13^23^ [—^^ ^ ) " (^34) 

Conversely, for weak wash-out in the r flavour, we have 



12 3MimiAmi f 390\ 

\Yb\<1T= 5 2-|^13|C23'\/ 777^ in^3 ^on ' (l^S) 



37g^ 32nv^ mf ' y mi V0.2 x 10^3 eV 589 

Thus, the maximal baryon asymmetry is obtained for intermediate wash-out effect. For 
definiteness we can estimate the upper bound on at — 3 x 10"'^ eV. For the smallest 
allowed value of mi for the QD spectrum, mi = 0.1 eV, we get: 

\Yb\ ^ 5 X 10-''\Ri,\^^^ . (136) 

Requiring that irfr <^ mi imposes a fine tuning on the values of the parameters Ru. 
For simplicity, we search for the solution of the equation frfr ~ 3 x 10"'^ eV ~ 0, which 
corresponds to: 

i?ll|f/rl| -i?12|f/r2|cOs(^)+i?13|f/.3|cOs(^) ^0 , (137) 

Ri2\Ur2\ Sin (^) - RislUrsl sin (^) ~ . (138) 

If a2i (tt3i) [032] = 0, we have that _Ri3 {R12) [Ru] = as well. Otherwise, the solution of 
Eqs. (11371) and (11381) is given by: 

fl„=ft3^5^!^ and R,, = R,,J^E^ (139) 

S12S23 sm a2i/2 C12S23 sm a2i/2 

Another possibility consists in having ttv — mi and small rrvz- However, also in this 
case the baryon asymmetry is suppressed due to the values of the CP-violating phases 
required to have <C mi, namely, a2i,a3i ~ 0,7r. In conclusion, for real R, it might be 
possible to reproduce the observed baryon asymmetry, but only for relatively large values 
of Ml ^ 10^^ GeV and small values of mi. A careful and detailed analysis should be 
performed on a case by case basis. A measurement of mi in the upper end of the range of 
allowed values [0.1 eV, 2.3 eV] in the quasi-degenerate spectrum would strongly disfavour, 
if not rule out, the possibility of having leptogenesis due uniquely to the low energy CP- 
violating phases, for hierarchical RH neutrinos and real R. 

For RiiRij = ±i\RiiRij\, the CP-asymmetry is enhanced by a factor 4mi/Am^. Also in 
this case an upper bound on |1b| can be derived, and it depends on rrv. A sufficiently large 
baryon asymmetry might be obtained for relatively large values of Mi. A more detailed 
analysis of this case is beyond the scope of the present study. 
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6 Baryon Asymmetry from Low Energy CP- Violating 
Dirac and Majorana Phases in U : the Case of Quasi- 
Degenerate RH Neutrinos 

In this Section we extend our previous findings to the case in which the RH neutrinos are 
quasi-degenerate in mass, Mi ~ M2 — M3 and the matrix R is real. We therefore consider 
the case in which the CP parities of the heavy and hght Majorana neutrinos are such that 
pf = p^j=l for all 2,j = 1,2,3. In such a case, indeed, CP invariance corresponds to 
having all the elements of the matrix R real, see Eq. ( l30i) . and 5 = 021 = 031 = 0. The 
degenerate pattern for the RH neutrino masses may arise if, for instance, there is a slightly 
broken SO (3) symmetry in the RH sector. The baryon asymmetry receives a contribution 
from the decay of all three RH neutrinos. The CP asymmetry in a given flavour / generated 
by the decay of the RH neutrino iV, (i = 1, 2, 3) is dominated by the one-loop self energy 
contribution [6H] and reads 



Sij = 72 — ' ^^l = - Ml , 

{AMl)' + MfT) ' 

4 = (AAt) .. (AAt) ^ {R^iR,irn,) ^ ^^K^,R,,R^M Wu) • (140) 

Notice, in particular, that = — Jjj and Sij = —Sji. The CP asymmetry e\ is resonantly 
enhanced when Tj = AMf^/Mi. At the resonance 



c Mj I 



(141) 



~ 2T.~ 2r, ' ^ 2 

The washing out of a given flavour / is now operated by the AL = 1 scatterings involving 
all three RH neutrinos. Therefore, the parameter is given by 

--E^-^-E-|C,/. (142) 

where we have set the nearly equal masses of the RH neutrinos approximately equal to M\. 
If the resonance operates for all three RH neutrinos, the CP asymmetry in the flavour I is 



- ^ Z^^ij - - 2^ hj - ^ 7— N / — ^ • (14^) 



i<j i<j (Ep ^pRip) (E, rngR% 
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It does not depend upon the mass of the RH neutrinos (if the running of the parameters 
is neglected). The elements of the matrix R can be parametrized by introducing a real 
antisymmetric matrix A 



R 



1 — cos r ^r. sin r , 
r H A . 
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^23 



^13 ) 



(144) 



where 1 is the 3x3 unity matrix and Aij = —Aji are the elements of the matrix A. In the 
limiting case where the RH neutrinos are exactly degenerate, one can perform an orthogonal 
rotation on the RH neutrino states which leaves the mass matrix of the RH neutrinos 
proportional to the unity matrix and defines a physically equivalent reparametrization of 
the RH neutrinos. This amounts to saying that for Mi = M2 = M3, the real matrix R can 
be set equal to the unity matrix (or A12 = ^23 = ^13 = 0). The flavour asymmetries in 
Eq. f ll40p vanish if i? = 1 which reflects the fact that no baryon asymmetry can be obtained 
in the exactly degenerate case. On the other hand, if the degeneracy is slightly broken, 
the elements of the matrix A are expected to be tiny, but not all vanishing and the baryon 
asymmetry is typically different from zero. We are now ready to study all three possible 
cases for the spectrum of the light neutrinos. We will restrict ourselves to the case in which 
the resonant condition is not satisfied for all couples of RH neutrinos and work under the 
condition that Mi ~ M2 ^ M3; the flavour asymmetries are generated resonantly only by 
the decays of A^i and N2 and we will set /13 ~ /23 — 0. 

We find 
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(145) 



in the normal hierarchical case; 
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in the inverted hierarchical case and, finally, 



X [{Ri2Rji - RiiR,2) Im iUaU:^) 

+ (i?i3i?i2 - i?i2i?j3) Im (fZ/af/fs) 

+ (i?,3/2ii - RiiRjz) Im (fZ/it/rs)] , 

mi c:^ m (147) 

in the degenerate case. In this latter case, since the washing-out factors are approximately 
same, the expressions for the baryon asymmetries may be simplified if R is real. Indeed, 
the total asymmetry ei vanishes and, if (10^ ^ Mi ^ 10^^) GeV, the fiavour asymmetry 
62 = — while rn2 — 2m ~ 2m^. One finds 

1 2 222 

Ys^- — —Y,. (148) 

In Fig. [T71 we show the correlation of the baryon asymmetry with the effective Majorana 
mass in neutrinoless double beta decay for the case of quasi-degenerate RH neutrinos 
and QD spectrum of light neutrinos. A number of projects aim to reach a sensitivity 
to \{m„)\ ~ (0.01 — 0.05) eV [69] and can certainly probe the region of values of 
for successfuU baryon asymmetry from the PMNS phases only. In particular, a direct 
information on the Majorana phase a2i may come from the measurement of (m^,), m, and 
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and might tell us if enough baryon asymmetry may be generated uniquely from the PMNS 
Majorana phases. 



7 Extension to the MSSM 

The extension to our findings to the supersymmetric version of the SM, the so-called Min- 
imal Supersymmetric Standard Model (MSSM) is rather straightforward [2S]. One has to 
consider the presence of the supersymmetric partners of the RH heavy neutrinos, the so- 
called sneutrinos Ni [i = 1,2,3), which also give a contribution to the fiavour asymmetries, 
and of the supersymmetric partners of the lepton doublets, the so-called slepton doublets. 
Since the effects of supersymmetry breaking may be safely neglected, the flavour CP asym- 
metries in the MSSM are twice those in the SM and double is also the possible channels by 
which a lepton flavour asymmetry is reproduced. However, the AL = 1 scatterings washing 
out the asymmetries are also doubled and the number of relativistic degrees of freedom is 
almost twice the one for the SM case. As a result, introducing new degrees of freedom 
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and interactions does not appreciably change the flavour asymmetries with respect to the 
values obtained within the SM. 

There are however two other and important differences with respect to the SM case. 
First, in the MSSM, the flavour-independent formulae can only be applied for tempera- 
tures larger than (1 + tan^/3) x 10^^ GeV, where tan/3 indicates the ratio of the vacuum 
expectation values of the two Higgs fields of the MSSM. Indeed, the squared charged lepton 
Yukawa couplings in the MSSM are multiplied by (1 + tan^/3). Consequently, charged /i 
and r lepton Yukawa couplings are in thermal equilibrium for (1 + tan^/3) x 10^ GeV <^ 
T -C (1 + tan^ f3) x 10^ GeV and all flavours in the Boltzmann equations are to be treated 
separately. For (1 + tan^ p) x 10^ GeV < T < (1 + tan^ p) x 10^^ GeV, only the r Yukawa 
coupling is in equilibrium and only the r flavour is treated separately in the Boltzmann 
equations, while the e and /i flavours are indistinguishable. This implies that the range of 
the RH (s)neutrino masses where flavour is relevant in leptogenesis is greater than the one 
in the SM by the factor (1 -|- tan^ j3) which is large even for moderate values of tan /5. As a 
consequence, the lower bounds given in Eqs. ( l93l) - ( l94l) change approximately to 

I sin ^13 sin sin > 0.11 {1 + tan^ p)-\ (150) 

|Jcp| ^ 2.4 X 10^^(1 +tan2/?)~^ (151) 

The lower bounds in Eqs. (l97l) - (l98l) and in Eqs. flllSp - flllQp change in a similar way. 
The shift of the range of the heavy Majorana neutrino masses, in which the lepton flavour 
effects are significant, to larger values has important implications also if the spectrum of 
the RH (s)neutrinos is hierarchical and the light neutrinos possess inverted hierachical (IH) 
spectrum. This is illustrated in Fig. [THl where we plot the baryon asymmetry versus the 
quantity Jcp for the IH spectrum of light neutrinos in the supersymmetric case for a given 
set of parameters. Let us recall that for real R matrix elements Ru and Ru, it was impos- 
sible to obtain baryon asymmetry compatible with the observations in the corresponding 
non-SUSY case. The reader should be also warned that, for values of tan f5 ^ 30, radiative 
corrections to the physical neutrino parameters should be accounted for |62j . 

Secondly, the relation between the baryon asymmetry and the lepton flavour asym- 
metries has to be modified to account for the presence of two Higgs fields. Between 
(1 + tan^ p) X 10^ and (1 + tan^ p) x 10^^ QeV, the relation is 



where the hat superscripts indicates that the flavour lepton asymmetries are computed 
including leptons and sleptons. Notice that if the spectrum of RH (s) neutrinos is quasi- 
degenerate as well as that of the light neutrinos, the wash-out factor are also the approx- 
imately same and the expressions for the baryon asymmetries may be simpliflied if R is 
real. Indeed, the total asymmetry ei vanishes and one of the flavour asymmetries may be 
expressed in terms of the others. Under these circumstances, one finds 

yMSSM ^ 10 447yMSSM QggN 

^ 31988 ^ ■ ^ ^ 
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Finally, in the case of supersymmetric leptogenesis one should also face the problem 
arising from the so-called gravitino bound. The latter is posed by the possible overpro- 
duction of gravitinos during the reheating stage after inflation. Being only gravitationally 
coupled to the SM particles, gravitinos may decay very late jeopardising the successfuU pre- 
dictions of Big Bang nucleosynthesis. This does not happen, however, if gravitinos are not 
efficiently generated during reheating, that is if the reheating temperature T^h is bounded 
from above, T^h < 10^° GeV [7D]. The severe bound on the reheating temperature makes 
the generation of the RH neutrinos problematic (for a complete study in the one-flavour 
case see [1]), if the latter are a few times heavier than the reheating temperature, rendering 
the thermal leptogenesis scenario difficult. There are, though, two possible ways out to this 
problem. First, leptogenesis might occur in a non-thermal way, that is the RH neutrinos 
might be generated not through thermal scatterings, but by other mechanisms, e.g. at the 
preheating stage [7T]. Alternatively, and maybe more interestingly, as we have previously 
seen if the two lightest RH neutrinos are quasi-degenerate in mass, the final baryon asym- 
metry does not depend upon their common mass. The latter, therefore, might be smaller 
than the largest possible reheating temperature and thermal leptogenesis might take place 
without any limitation from the gravitino bound. 



8 Conclusions 

In this paper we have systematically investigated the connection between the leptogenesis 
and the low energy CP- violation in the lepton (neutrino) sector. Our study was stimulated 
by the recent progress in the understanding of the importance of lepton flavour effects in 
leptogenesis. It lead to the realization that these effects can play a crucial role in the 
leptogenesis scenario, both from the quantitative and the qualitative point of view. When 
the lepton flavour effects are taken into account, the flnal baryon asymmetry is the sum 
of three different contributions given by the CP asymmetries generated in each flavour 
(lepton charge), properly weighted by the corresponding wash-out factor. In the one- 
flavour approximation, which holds only if leptogenesis is taking place at a temperature 
higher than about 10^^ GeV, the final baryon asymmetry is proportional to the total baryon 
CP asymmetry (summed over the three flavours) and weighted by a single wash-out factor 
(obtained by summing the wash-out factor of the three lepton flavours). 

There are many differences between the predictions for the baryon asymmetry Yb ob- 
tained in the one-flavour approximation and in the case when the flavour effects are ac- 
counted for. The baryon asymmetry Yq, derived in the one-flavour approximation, for 
instance, vanishes if the light neutrinos are degenerate in mass. Correspondingly, Yq has to 
be proportional to a difference of masses of the light neutrinos. In the "flavour" case this 
suppression can be absent even when the leptogenesis CP violation is due entirely to the 
low energy phases in the PMNS matrix. However, the most signiflcant difference is that 
in the one flavour approximation there is no direct connection between the leptogenesis 
CP-violating parameters and the CP-violating parameters - Dirac and Majorana phases, 
present in the lepton (neutrino) sector. In particular, a possible future observation of CP 
violation in neutrino oscillations would not automatically imply, within the "one-flavour" 
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leptogenesis scenario, the existence of a baryon asymmetry. In the "flavoured" treatment 
of leptogenesis, however, this conclusion does not universally hold and the observation of 
CP violation in the lepton (neutrino) sector would generically imply a nonvanishing baryon 
asymmetry. Including the effects of lepton flavour, therefore, allows to build a new bridge 
between the CP violation in leptogenesis and the observables depending on the CP-violating 
Dirac and Majorana phases in the PMNS neutrino mixing matrix, such as the CP violating 
rephasing invariant Jqp which controls the magnitude of CP-violation effects in neutrino 
oscillations, the effective Majorana mass \{m) \ in neutrinoless double beta decay, etc. The 
study of such a connection has been the main subject of our paper. 

We have first derived the constraints the requirement of CP-invariance imposes on 
the neutrino Yukawa couplings A and on the elements of the complex orthogonal matrix R 
appearing in the "orthogonal" parametrisation of A. The CP-parities of the light and heavy 
Majorana neutrinos, which take the values ±z, play a special role in these constraints. The 
CP-invariance constraints are useful for understanding the source of CP violation generating 
the CP asymmetries in the heavy Majorana neutrino decays. One example is the case of real 
matrix R and specific CP-conserving values of the Majorana and Dirac phases in the PMNS 
neutrino mixing matrix, which corresponds to violation of CP-symmetry at high energy 
in leptogenesis, leading to the generation of non-zero baryon asymmetry. The indicated 
constraints help to clarify also under which conditions the leptogenesis CP-asymmetries 
are due entirely to the low energy CP-violating phases of the PMNS matrix. 

Taking into account the lepton flavour effects in leptogenesis, we have subsequently 
investigated in detail the possibility that the CP-violation necessary for the generation of 
the baryon asymmetry of the Universe is due exclusively to the Dirac and/or Majorana 
CP-violating phases in the PMNS matrix, and thus is directly related to the low energy 
CP- violation in the lepton sector (e.g., in neutrino oscillations, etc.). We have derived 
results for two types of spectrum of the heavy RH Majorana neutrinos: i) hierarchical, in 
which the lightest RH neutrino A'^i is much lighter than the other two RH neutrinos, and 
ii) quasi-degenerate, in which the two lightest RH neutrinos Ni^2 are almost degenerate in 
mass and have masses which are smaller than the mass of the third one. For each of the two 
cases, we have presented predictions for the baryon asymmetry for three types of spectra 
of the light Majorana neutrinos: normal hierarchical (NH), inverted hierarchical (IH) and 
quasi-degenerate (QD). In all numerical calculations we have used the best fit values of the 
solar and atmospheric neutrino oscillation parameters, Am\, Am^, sin^ ^23 and sin^6'i2, 
given in Section 2. 

For hierarchical RH neutrino mass spectrum, the lepton flavour effects are relevant in 
leptogenesis for Mi ^ xlO^^ GeV. The predicted baryon asymmetry 1b depends linearly 
on Ml. In order to reproduce the observed baryon asymmetry, generically values of Mi ^ 
3 X 10^° GeV are required. We have shown that if the light neutrinos have a NH spectrum, 
the requisite baryon asymmetry can be produced if the only source of CP violation is either 
the Majorana phases or the Dirac phase in the PMNS matrix, ?7pmns (Figs- [D-E]). When the 
only CP- violating parameter is the low energy Majorana phase 032 = 031 —a2i) we can have 
successful leptogenesis as long as | sin(a;32/2)| is not exceedingly small and Mi ^ 3.5 x 10^° 
GeV. If the only source of CP violation is the Dirac phase 6 in f/pMNS, the observed baryon 



44 



asymmetry can be reproduced provided Mi ^ 2 x 10^^ GeV and | sin^ia sin(5| ^ 0.1, ^13 
being the CHOOZ angle. This condition leads to the inequality sin 6*13 ^ 0.1 and to the 
following lower bound on the CP-violating rephasing invariant Jcp, associated with the 
Dirac phase in the PMNS matrix: |Jcp| ~ 2.0 x 10"^. Values of sin^is ^ 0.1 can be 
probed in the forthcoming Double CHOOZ and Daya Bay reactor neutrino experiments. 
CP- violation effects with magnitude determined by |Jcp| ~ 2.0 x 10"^ are within the 
sensitivity of the next generation of neutrino oscillation experiments, designed to search 
for CP- or T- symmetry violation in the oscillations. Moreover, since in this case both 
IIbI oc I sin ^13 sin 5 1 and |Jcp| | sin ^13 sin 5 |, given the other parameters on which |1b| 
and I Jcp I depend, there exists a correlation between the rephasing invariant Jcp, which 
controls the magnitude of the CP- violation effects in neutrino oscillations, and the baryon 
asymmetry (Fig. [U]). 

In the case of IH light neutrino mass spectrum and negligible lightest neutrino mass 
777,3, the observed baryon asymmetry cannot be reproduced if the product of the elements 
of the matrix i?, RuRu, is purely real: the generated baryon asymmetry is generically 
small, being suppressed by the additional factor Amg/Am^. However, if R11R12 is purely 
imaginary (and CP-conserving), a sufficiently large baryon asymmetry compatible with the 
observations can be obtained both when the only source of CP-violation is the Majorana 
phase a2i, or the Dirac phase 6, in ?7pmns (Figs- [7]-[TnD. In the case of Majorana CP- 
violation, depending on the sgn(Im(i?iii?i2)), values of Mi ^ 5 x 10^° GeV or somewhat 
larger (e.g.. Mi ^ 1.6 x 10^^ GeV) are required. Since both the baryon asymmetry \Yb\ and 
the effective Majorana mass in (/3/3)oj/-decay, |(m)|, depend on the Majorana phase 021, for 
given values of the other parameters there exists a direct correlation between the values 
of \Yb\ and \{m)\ (Fig. [TU]). We have shown that one can have successful leptogenesis in 
the case under discussion also if S13 7^ and the CP- violation is generated only by the 
Dirac phase 6 in ?7pmns (Figs. fTTl and fT2i) . For Mi ^ 5 x 10^^ GeV, the observed baryon 
asymmetry can be reproduced if | sin^i3 sin(5| ^ 0.02. This requirement implies that we 
should have also sin ^13 ^ 0.02 and |Jcp| ~ 4.6 x 10~^. Values of sin ^13 and of |Jcp| as 
small as 0.02 and 4.6 x 10~^, respectively, can be probed in neutrino oscillation experiments 
at neutrino factories. There exists a correlation between the rephasing invariant Jcp and 
the baryon asymmetry Yb in this case as well (Fig. [T3l) . 

The analysis we have performed showed that if the light neutrinos have QD spectrum, 
the baryon asymmetry is generically too small mainly due to the large wash-out suppression 
factor. 

For heavy RH Majorana neutrinos with QD spectrum, leptogenesis takes place through a 
resonance effect. The main new feature is that the final baryon asymmetry does not depend 
on the mass of the RH neutrinos. This property is crucial, allowing the generation of a 
sufficiently large baryon asymmetry even in the case of QD light neutrino mass spectrum. 
In the latter case the predicted baryon asymmetry is correlated with the effective Majorana 
mass in the neutrinoless double beta decay (Fig. [T71) . 

Finally, we have discussed how the results on leptogenesis we have obtained will be 
modified in the minimal supersymmetric extension of the Standard Theory (MSSM) with 
right-handed Majorana neutrinos and see-saw mechanism of neutrino mass generation. We 
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have noticed, in particular, that for hierarchical heavy Majorana neutrino mass spectrum, 
the range of the lightest RH neutrino mass Mi for which the lepton flavour effects are 
relevant in leptogenesis is greater than the one in the non-supersymmetric case by the factor 
(1 + tan^ P), Ml ^ (l+tan^ P) x 10^^ GeV, tan j3 being the ratio of the vacuum expectation 
values of the two Higgs fields of the MSSM. This can have important implications especially 
in the cases when the generation of the baryon asymmetry in the non-supersymmetric case 
is strongly suppressed. We have also stressed that a quasi-degenerate spectrum of the 
heavy RH Majorana neutrinos is welcome in the case of supersymmetric leptogenesis since 
it renders the gravitino bound harmless. 

The results obtained in the present article underline the importance of understanding 
the status of the CP-symmetry in the lepton sector and, correspondingly, of the experiments 
aiming to measure the CHOOZ angle 6*13 and of the experimental searches for Dirac and/or 
Majorana leptonic CP-violation at low energies. 
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Figure 1: The baryon asymmetry Yg as a function of R12 in the case of real R12 and -R13, 
sign (_Ri2-Ri3) = +1 (/323 = 0), + -R13 = 1, ■§13 = 0.20, hierarchical RH neutrinos and 
NH light neutrino mass spectrum and a) Majorana CP- violation (blue line), 5 = and 
«32 = 7r/2 (k = +1), and b) Dirac CP-violation (red line), 5 = 7r/2 and 032 = (k' = +1), 
for Ml = 5 X 10^^ GeV. The neutrino oscillation parameters Amg, sin^^i2, Am\ and 
sin^ 2^23 are fixed at their best fit values. 
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Figure 2: The same as in Fig. [T]but for sign (_Ri2-Ri3) = —1 {P23 = tt) and a) Majorana 
CP- violation (blue line), 6 = and 032 = 7r/2 (k = —1), and b) Dirac CP- violation (red 
line), 6 = 7r/2 and 032 = (k' = —1). 
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Figure 3: The baryon asymmetry as a function of the Majorana phase varying in the 
interval 0:32 = [0, 27r] in the case of Majorana CP-violation, hierarchical RH neutrinos and 
NH light neutrino mass spectrum, for 6 = 0, real R12 and R13, \Ri2\ = 0.92, |i?i3| = 0.39, 
sgn(i?i2-Ri3) = +1 (/?23 = 0, K = +1), Ml = 5 X 10^° GeV, and two values of S13: S13 = 
(blue line) and 0.2 (red line). 
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Figure 4: The same as in Fig. [3] but for real R12 and R13 having opposite signes, 
sgn(i?i2-Ri3) = —1 {P23 = n, K = —1), |i?i2| = 0.92, |i?i3| = 0.39, and two values of 
Si3- Sis = (red line) and S13 = 0.1 (blue line). 
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Figure 5: The baryon asymmetry \Yb\ as a function of the Dirac phase 6 varying in the 
interval 6 = [0, 27i] in the case of Dirac CP- violation, 032 = 0; 27r, hierarchical RH neutrinos 
and NH light neutrino mass spectrum, for Mi = 5 x 10^^ GeV, real Ru and -R13 satisfying 
\Ri2\^ + \Ri3\^ = 1, 1^12! = 0.86, |i?i3| = 0.51, sign (i?i2i?i3) = +1, and for i) 032 = 
(k' = +1), S13 = 0.2 (red line) and S13 = 0.1 (dark blue line), ii) 032 = 27r (k' = — 1), 
S13 = 0.2 (light blue line). 
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Figure 6: The correlation between the rephasing invariant Jcp (in blue) and the baryon 
asymmetry Yb when varying the Dirac phase 6 = [0,27r], in the case of hierarchical RH 
neutrinos and NH light neutrino mass spectrum and for S13 = 0.2, 032 = (27r), I-R12I = 
0.86, |i?i3| = 0.51, sign (i?i2i?i3) = +1 (-1) i/323 = (tt), k' = +1), Mi = 5 x lO^^ GeV . 
The red region denotes the 2a allowed range of 1b- 
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Figure 7: The baryon asymmetry Yb as a function of |-Rii| in the case of hierarchical 
RH neutrinos and IH light neutrino mass spectrum, Majorana CP-violation, 6 = and 
Q^32 = 7r/2, Ml = 2 X 10^^ GeV, purely imaginary R11R12 = inlRuRul and k, = +1 (dark 
blue and red lines), k = — 1 (light blue and green lines), |i?i2 
(green and red lines) and S13 = (light and dark blue lines). 



^i2p-|i?i3p = 1, andfor si3 = 0.2 
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Figure 8: The baryon asymmetry as a function of the Majorana phase 032 varying in the 
interval = [0, 27r] in the case of hierarchical RH neutrinos and IH light neutrino mass 
spectrum, Majorana CP-violation, 6 = 0, purely imaginary RuRu = ifi:|-Rii-Ri2|, n = +1 
{P12 = 7r/2), - |i?i2p = 1, |i?ii| = 1.05, and for Mi = 2 x 10" GeV, and two values 

of S13: S13 = (blue line) and 0.2 (red line). 
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Figure 9: The same as in Fig. [HI but for k = — 1 (/5i2 = 37r/2) and |-Rii| = 1.2. 
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Figure 10: The baryon asymmetry \Yb\ versus the effective Majorana mass in neutrinoless 
double beta decay, in the case of Majorana CP- violation, hierarchical RH neutrinos 

and IH light neutrino mass spectrum, for 6 = 0, S13 = 0, purely imaginary R11R12 = 
iK\RiiRi2\, K = +l (/5i2 = 7r/2), ji^np - {R-^^l^ = 1, \Ru\ = 1.05 and Mi = 2 x 10^^ GeV. 
The Majorana phase 021 is varied in the interval [— 7r/2, 7r/2]. 




^11 

Figure 11: The baryon asymmetry \Yb\ as a function of |-Rii| in the case of hierarchical RH 
neutrinos and IH light neutrino mass spectrum and Dirac CP- violation, 6 = it/2, a2i = vr 
(k' = -|-1), purely imaginary R11R12 = iK|-Rii-Ri2| (I-R12P — l-RisP = 1), for k = — 1 (red 
and dark blue lines) and k = +1 (light blue and green lines), S13 = 0.2 (light blue and dark 
blue lines) and S13 = 0.1 (green and red lines), and Mi = 2 x 10^^ GeV. 
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Figure 12: The asymmetry \Yb \ as a function of the Dirac phase 6 in the case of hierarchical 
RH neutrinos, IH light neutrino mass spectrum, Dirac CP-violation, a2i = tt [k' = +1), 
RiiRu = i K |i?ii_Ri2| (l-RiiP — I-R12P = 1)) ^ = ^1 (red and dark blue lines), k = +1 
(light blue and green hues), for Mi = 2 x 10^^ GeV, and S13 = 0.1 (red and green hues) 
and Si3 = 0.2 (dark blue and light blue lines). Values of |-Rii|, which maximise \Yb\ have 
been used: |-Rii| = 1.05 in the case of k = —1, and |-Rii| = 1.3 (1.6) for k, = +1 and 
si3 = 0.2 (0.1). 
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Figure 13: The correlation between the rephasing invariant Jqp (in blue) and the asym- 
metry Yq in the case of hierarchical RH neutrinos, IH hght neutrino mass spectrum, Dirac 
CP-violation, 021 = vr (k' = +1), RnRu = i k |i?iii?i2| (l-RiiP — I-R12P = 1), = +1, 
and for S13 = 0.2, Mi = 5 x 10^° GeV and |-Rii| = 1.3. The Dirac phase 5 is varied in the 
interval [0, 27r]. The red region denotes the 2cr allowed range of 1b- 



57 







1 1.2 1.4 1.6 1.8 2 2.2 

^11 



Figure 14: The baryon asymmetry Yb as a function of |-Rii| in the case of hierarchical RH 
neutrinos and IH hght neutrino mass spectrum, 6 = 7r/2, 021 = 0, and purely imaginary 
-R11-R12 = ifi\RiiRi2\ satisfying I-R12P — l-RisP = 1, for k = +1 (red and dark blue lines), 
K = —1 (light blue and green lines). Mi = 2 x 10^^ GeV and two values of S13: S13 = 0.2 
(red and light blue lines) and S13 = 0.1 (dark blue and green lines). 




Figure 15: The baryon asymmetry as a function of the Dirac phase 6 in the case of 
hierarchical RH neutrinos and IH light neutrino mass spectrum, for 021 = 0, Mi = 2 x 10^^ 
GeV, purely imaginary R11R12 = i n |i?iii?i2| (l-RiiP — I-R12P = 1); with k = +1 (green 
and light blue lines ) and k = —1 (red and dark blue lines), and for |-Rii| = 1.05 and two 
values of S13: S13 = 0.1 (light and dark blue lines), S13 = 0.2 (red and green lines). 
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Figure 16: The correlation between the rephasing invariant Jcp (in blue) and the baryon 
asymmetry Yb in the case of hierarchical RH neutrinos and IH light neutrino mass spectrum 
for = 0, si3 = 0.2, Ml = 2 X 10" GeV, RuRu = i {RiiRul i\Rn\^ - |^i2p = 1), 
and for k = +1 and |-Rii| = 1.05. The Dirac phase S is varied in the interval [0,27r]. The 
red region denotes the 2a allowed range of Yq. 
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Figure 17: The baryon asymmetry versus the (/3/9)oi/-decay effective Majorana mass 
in the case of QD heavy RH neutrinos and QD light neutrino mass spectrum, and for 
6 = n/3, Si3 = 0.01, Ml = 10^^ GeV and m = 0.1 eV. The Majorana phase 032 is varied in 
the interval [0,7r/3]. 
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Figure 18: The rephasing invariant Jqp versus the baryon asymmetry in the case of 
super symmetric hierarchical RH neutrinos and IH hght neutrino mass spectrum and for 
= 7r/4, si3 = 0.1, Ri2 = 0.86, R13 = 0.5, sign {RuRis) = +1 and Mi = 6 x 10^^ Q^y 
The Dirac phase 6 is varied in the interval [0,27r]. 
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